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ABSTRACT 

> 

cn 

' In this paper we present mid- infrared spectra of a comprehensive set of Herbig Ae/Be 

■ stars observed with the Spitzer Space Telescope. The signal-to-noise ratio of these 

QQ . spectra is very high, ranging between about a hundred and several hundreds. During 

O , the analysis of these data we tested the validity of standard protoplanetary dust models 

and studied grain growth and crystal formation. On the basis of the analyzed spectra, 
the major constituents of protoplanetary dust around Herbig Ae/Be stars are amorphous 
silicates with olivine and pyroxene stoichiometry, crystalline forsterite and enstatite and 
• silica. No other solid state features, indicating other abundant dust species, are present 

^ in the Spitzer spectra. Deviations of the synthetic spectra from the observations are 

most likely related to grain shape effects and uncertainties in the iron content of the 
dust grains. 

Our analysis revealed that larger grains are more abundant in the disk atmosphere of 
flatter disks than in that of flared disks, indicating that grain growth and sedimentation 
decrease the disk flaring. We did not find, however, correlations between the value of 
crystallinity and any of the investigated system parameters. Our analysis shows that 
enstatite is more concentrated toward the warm inner disk than forsterite, in contrast 
to predictions of equilibrium condensation models. None of the three crystal formation 
mechanisms proposed so far can alone explain all our findings. It is very likely that all 
three play at least some role in the formation of crystalline silicates. 
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Subject headings: circumstellar matter - infrared:planetary systems - infrared:stars - 
stars:formation - stars:pre-main-sequence 



1. Introduction 



The class of Herbig Ae/Be (hereafter HAeBe) stars was estabhshed by Herbig (jl96ol ) as stars 
which are surrounded by nebulosities and the optical spectra of which show emission lines. Further 
investigations revealed that these sources are young stars (1-10 Myr) with masses between 2 and 
10 Mq in the later stages of their pre-main sequence evolution. Observations of these sources 
at infrared wavelengths revealed that excess er nission above the stellar p hotosphere is another 
characteristic of HAeBe stars (for a review, see IWaters &: Waelkena (119981 )). The infrared excess 
emission arises from a protoplanetary disk (jWaters fc Waelkena Il998l ) and in many cases from 



an envelope as well (jLeinert et al 



2001 



HAeBe stars are, therefore, frequently regarded as the 
higher mass counterparts of the low- mass T Tauri stars. Planet formation theories suggest that 
this evolutionary stage (l~10Myr), is exactly where the formation of planetary embryos is likely 
to occur. Thus, HAeBe stars are natural candidates for studying the physical processes playing an 
important role in planet formation. Stars of spectral type A have gained a rene wed interest becaus e 



of the recent direct imaging detection of extrasolar planets around these stars (iMarois et al.ll2009l ) 



In this study, we focus on a subgroup of the HAeBe class with spectral type between late 
B and A-F, i.e., the lower mass end of the HAeBe class(hereafter HAe stars). Spectral energy 
distributions (SEDs) of these sta rs can be well represent ed with models of a passive protoplanetary 



disk with a puffed- up inn er rim (IDullemond et al.ll200ll ). Based on observations with the Infrared 



Space Observatory (ISO). iMeeus et al.l (120011 ) classified the HAe stars into two groups. SEDs of 
Group II sources can be well fitted with a power law at mid- to far-infrared wavelengths. An 
additional blackbody component is required, however, to fit the SEDs of Group I sources at far- 
infrared wavelengths. Theoretical models of protoplanetary disks showed that SEDs of Group I 
sources can be explained by flared disks, which are in vertical hydrostatic equilibrium and where 
gas and dust are well mixed. Later on as dust grains grow in size and settle to the mid-plane, the 
disk becomes flatter producin g the steeper, bluer mid- to far- infrared SEDs of Group II sources 
(jPullemond Dominikll2004al ). 



The global shape of the SED, however, carries only limited information about the physical prop- 
erties of protoplanetary dust grains and the processes they undergo. Mid-infrared spectroscopy, on 
the other hand, is an excellent diagnostic tool for studying the size, shape and chemical composition 
of protoplanetary dust grains. The mid-infrared domain is rich in vibrational resonances of silicates 
with different compositions, which are the main constituents of the protoplanetary dust (Henning 
2009). These mid- infrared emission features originate in the hot surface layer of the disks where 
the temperature is above ~100K. Since this region of the disk is optically thin, by analyzing the 
emergent spectra, the composition of the dust mixture as well as the physical parameters of the 
radiating dust grains (e.g., size or shape) can be derived. Mid-infrared spectroscopy, however, has 
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also limitations. It is sensitive only to dust grains which show resonances in the mid-infrared, i.e., 
grains larger than several microns or "featureless" grains (e.g., amorphous carbon or iron) cannot 
be studied in this way. Since mid-infrared features arise from the surface layers of the disk, the 
derived grain properties are not necessary representative for the whole vertical extent of the disk. 



Crystalline silicates are abundant in many solar system comets (see e.g.. lWooden et al.l (120071 ) 
and references therein) but they are esse ntially missing f rom t he int erstellar m e dium ( ISM). Prom 
the analysis of the 10 silicate band iKemper et al.l (|2005l ) and iMin et al.l (|2007al ) placed an 
upper limit of 2 % in terms of mass for the abundance of silicate crystals in the ISM. Although they 
represent usually a minor dust constituent in terms of abundance, the sharp featu res of crystalline 



silicat es are frequently observe d toward young stars, including HAe stars (e.g. iBouwman et al 



( 2001 ). van Boekel et al. (j2005l )). It is therefore reasonable to assume that crystallization occurs in 



the disks of young stars. Due to their sharp features crystalline silicates can be used as tracers to 
investigate dynamic processes in protoplanetary disks. Both ways of crystal formation (annealing 
and direct condens ation from the gas phase) require high temperature, typically above 1000 K 
(jFabian et al.ll2000l ). The fact that we still observe crystals in the outer disk where the temperature 
is of t he order of 100 K suggests, that either a large-scale mixing should occur in the disks of young 
stars (lGailll2004l) or amorphou s grain s should be hea t ed loca lly (e.g. shocks) to be transformed into 
crystals (e.g. iHarker fc DeschI (|2002l ^. ISargent et al.l (|2009al )^. 



Grain growth is another important process in protoplanetary disks which can be studied by 
mid-infrared spectroscopy. As sub-micron sized amorphous grains grow in size above a micron 
their 10 fim silicate feature becomes broader and flatter compared to the triangular shaped feature 
of the smaller grain population. Thi s was indeed observed i n the spectrum of m any young stars 



regard l ess of their spe c tral t ype (e.g. Ivan Boekel et al.l (120051) : I Apai et al.l (120051 ): iBouwman et al 



Bouwman et al 



(2003) and 



Meeus et al 



(2008): IWatson et al.l (|2009l )). It was also reported by 
(120091 ) that the size of the dust grains tends to be larger in flatter disks compared to flared one. 
This is the first observational evidence that dust sedimentation can be the reason why initially 
flared disks evolve to flatter ones. 

In all of the above mentioned studies the average signal-to- noise ratio (S/N) of the mid- infrared 
spectra was of the order of ~100 or lower. In this paper, we take one step further and analyze 
Spitzer IRS spectra of a comprehensive set of HAe stars with extremely high quality (with S/N up 
to several hundreds). The goal of our analysis is to (1) test our knowledge collected from analysis 
of lower S/N data and (2) look for possible new dust species/effects, which are rare/weak enough 
to be observable only in high quality da t a, (3) i nvestigate if the re latioi iship between disk fla ring 
and grain size, found by lBouwman et al.l (j2008l ). iMeeus et al.l (j2009l ) and ISarg ent et al.l (l2009bi ) for 
T Tauri stars, also holds for HAe stars. The analysis of the PAH emission will be presented in a 
separate paper (Acke et al., in prep). 
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Observations 



2.1. Sample selection 



The list of source s was c ompiled from the samples oflThe et alj (Il994l). ISvlvester et al.l (| 19961 ) . 
van den Ancker et al.1 (|l998l ). ISvlvester fc Manningsl (|20od ) and lMalfait et alJ (|l998al ). This source 
list was cross- ;^.d— h the Spitzer Archival andlTob.ervlI sources were selected. Since the 
five studies, from which our original sample was derived, used different classification criteria HAeBe 
stars, the sample was not uniform and false classifications occurred. Sources which were in fact not 
Herbig Ae stars (but e.g., classical Be systems or asymptotic giant branch stars) were rejected from 
our sample. Sources with obvious extended emission were also rejected, since our goal was to study 
the dust properties in disks around HAe stars. The resulting final sample consists of 53 sources 
in total out of which 45 shows silicate emission features while the remaining eight sources show 
only emission lines from polycyclic aromatic hydrocarbons (PAHs). Although in some "PAH-only" 
sources there may be a hint of a weak 10 /im silicate feature the feature is so weak/shallow that 
a meaningful dust composition cannot be determined we therefore excluded them from the final 
sample. The final list of targets is presented in Table[TJ In this paper we focus on the sources 
showing silicate emission and the analysis of the PAH bands will be presented in a forthcoming 
paper (Acke et al., in prep). 



2.2. Data reduction 



The spectra presented in this paper were obtained using the Infrared Spectrograph (IRS lHouck et al 



20041 ) on board the Spitzer Space Telescope. In most cases observations were performed using the 
short-low module (5.2-14.5 /im) of the low-resolution (R=~60-120) spectrograph and both the 
short-high (9.9-19.5 fim) and long-high (18.7-37.2 fim) modules of the high resolution (R=600) 
spectrograph. In the case of HD152404 only low-resolution modules, short-low and long-low (14- 
35 //m), were used. For 8 sources there were no low-resolution spectra taken with the Spitzer IRS 
instrument, only the short-high and the long-high modules were used. 

In the case of low-resolution mode the data reduction process started from the droopres 
intermediate data product processed through the SSC pipeline S15.3.0. Our data are further 
processed using spectral extraction tools devel oped for the PEPS Spitzer science legacy program, 



partially based on the SMART software package (jHigdon et al.ll2004l ). Most of our observations were 



taken in standard staring mode where the target is observed at the two nominal nod positions in 
the slit (~18" from the slit center), using multiple cycles per target for redundancy and to allow 
the rejection of artifacts introduced by bad pixels or cosmic ray hits. A high accuracy IRS or 
PCRS peak-up (with a la pointing uncertainty of 0.4" radius) was used to acquire targets in the 



^Most of the observations came from two programs, PI: J. Bouwman, PID:3470 and PI: B. Acke, PID:20308 
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spectrograph slit. A subset of our sources has been observed in 2x3 mapping mode without a 
peak- up. The smah maps consist of two positions at the nominal nod positions in slit, similar 
to normal staring observations, and three map positions in a perpendicular direction to the slit, 
with the central position centered on the target and the other positions shifted by half a slit width 
(1.8"). Effectively, this results in three standard staring mode observations with one observation 
reasonably centered on the source and two offset observations. We have used the central map 
position and use those as normal nodded observations in standard staring mode. 

As a first step, we correct for the background emission and stray light (mainly coming from 
the infrared background seen by the peak-up array) by subtracting the associated pairs of imaged 
spectra of the two nodded positions along the slit for each module and order. Pixels flagged by the 
data pipeline as being "bad" were replaced with the average pixel value of a six pixel elongated 
box surrounding the bad pixel. The method we apply for finding the mean pixel value resembles 



Nagao & Matsuyama filtering (iNagao &: Matsuvamal 11973 ) and ensured edge preservation in the 



source region of our spectral images. The spectra were extracted using a fixed-width aperture of 
six pixels centered on the position of the source. The exact source position relative to the slit was 
determined by fitting a sine profile to the spectra in the dispersion direction using the collapsed 
and normalized source profile. 

The spectra are calibrated with a relative spectral response function derived from IRS spectra 
and MARCS stellar models for a suite of calibrators provided by the Spitzer Science Center through 
the Spitzer data archive. The spectra of the calibration stars (ryl Dor, HR 6606, HR 7341) were 
extracted using the same method as for our science targets. One of the most difficult problems with 
spectroscopy using a narrow slit is the spectro-photometric calibration. Due to telescope pointing 
uncertainties and drifts, a variable fraction of source flux is being blocked by the slit. For high 
accuracy peak-up observations the intrinsic photometric accuracy is about 10%, while observations 
with no peak-up have a far lower accuracy. Due to the wavelength dependence of the point spread 
function (PSF) these pointing-induced flux losses will also change the spectral shape. To remove 
any effect of pointing offsets, we developed a correction method based on the PS F of the IRS 



i nstru ment, correcting for possible flux losses. For details of this method we refer to ISwain et al 



(|2008l ). We estimate the flux accuracy we can achieve with our data using this method to be 



The data reduction procedure for the hi gh-resolution data w as based on the method developed 



by the Cores-to-Disks Spitzer legacy team (|Lahuis et al.l 120071 ) . The procedure started from the 



rsc products processed through the same version (S15.3.0) of the Spitzer data pipeline as the low- 
resolution data. The spectra were extracted in two ways. The first method uses a fixed width 
aperture very similar to the method we used for the low-resolution data. The second method 
is an optimal source profile extraction method which fits an analytical PSF derived from sky- 
corrected calibrator data and an extended emission component, derived from the cross-dispersion 
profiles of the flat-field images, to the cross-dispersed source profile. It is not possible to correct 
for the sky contribution in the high-resolution spectra, subtracting the two nod positions as with 
the low-resolution observation, due to the small slit length. We either subtracted an observation 
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on the sky at a position close to the source or, when no such sky observation was taken, used 
the background estimate from the source profile fitting extraction method. For correcting "bad" 
pixels we used the I RSCLEAN package. We further removed low-level (~1%) fringing using the 



irsfringe package (jLahuis Sz Boogertl l2003l ) . We carefully checked that our fringe removal was 
not affecting the multiple silicate bands seen in our spectra. As the frequency of the fringes is 
reasonably well constrained and higher than the typical width of the observed thermal emission 
features from the various dust components, we found this not to be a problem. 

The flux calibration for the high-resolution spectrograph has been done in a similar way as for 
the low-resolution observations. For the relative spectral response function we also used MARCS 
stellar models and calibrator stars provided through the Spitzer Science Center. The spectra of the 
calibration stars were extracted in an identical way to our science observations using both extraction 
methods. As with the low-resolution observations, we also corrected for possible flux losses due 
to pointing offsets. We estimate the absolute flux calibration uncertainty for the high resolution 
spectra to be ~3%, slightly higher than that of the low-resolution observations. We found that the 
fixed width aperture extraction gave the best result for the short-high module (9.9-19.5 /^m), while 
for the long-high module (18.7-37.2 /um) the optimal extraction method was slightly better. For 
the final spectra presented in this paper we therefore used the results of the fixed width aperture 
extraction for the wavelengths shortward of 19 /im, and the results of the optimal extraction method 
for wavelengths longward of 19 /im. We want to note that in a few spectra (HD35187, HD38120, 
HD139614) the spectrum in the 12th order (~31-34 ^m) of the long-high module seems to be tilted 
compared to the neighboring orders using the optimal extraction method. In the case of HD139614 
we saw similar behaviour in the 15th order of the long- high module (~25-27 ^m). Although the full 
aperture extraction method did not show such strong tilt, it gave significantly higher noise level in 
this order, than the optimal extraction. Since the choice of the extraction method did not change 
our results we used the optimal extraction method for the long-high module to obtain uniformly 
reduced data in the whole sample. 

After the spectra have been reduced the different modules were combined to achieve our final 
spectra. Between 5.5 and 13.5 /um the short-low module was used while we used the short-high and 
the long-high for the 13.5-19.5/im and 19.5-35 micron wavelength intervals, respectively. For the 
sources, where no low-resolution Spitzer IRS spectra were taken, Spitzer spectra were supplemen ted 



shortward of 13.0 //m by data taken with the TIMMI2 instrument from Ivan Boekel et al. 



(12005!), if 



such data were available. The high-resolution and the TIMMI2 spectra were rebinned for a uniform 
spectral resolution of R=160 for the spectral fitting. Though the absolute flux calibration of the 
IRS observations is very good, any differences in the absolute flux calibration in various modules 
were handled in the following way. The spectra in different modules were scaled to a reference 
module which is chosen to be the one with the lowest absolute flux calibration uncertainty. We 
used, therefore, the short-low module as a reference, if it was present. If no short-low module was 
available the short-high module was chosen to be the reference. The applied scaling factors are of 
the order 1. 
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3. Analysis 
3.1. Dust model 

In order to study evolution and thermal processing of protoplanetary dust grains using mid- 
infrared spectroscopy, first one needs to identify the abundant dust species in the disks around 
young stars. Such an identification can be done by comparing the laboratory measurements of 
mass absorption coefficients (MACs) of different materials to the emission fe atures observe d in th e 
spectra. Such a comparison/identification has already been done by e.g., iMolster et al.l (|2002l ). 
These studies showed that mid-infrared spectra of young stars can be well reproduced by a mix- 
ture of five dust species, amorphous silicates with olivine and pyroxene stoichiometry, crystalline 
forsterite and enstatite and silica. The IRS instrument on board the Spitzer Space Telescope al- 
lowed us to improve the S/N of mid-infrared spectra by more than an order of magnitude compared 
to ISO SWS and and by a factor of 5-8 compared to ground-based instruments (e.g., COMICS, 
TIMMI2, T-ReCS). The exercise was repeated on the Spitzer data and emission features seen in 
the spectra were identified. The identification of the features is summarized in Tab[3l The dust 
features seen in our spectra can be identified as any of the following materials: amorphous silicates 
with olivine and pyroxene stoichiometry, forsterite, enstatite, silica and PAHs. 

Amorphous silicates of olivine (MgxFei_xSi04) and pyroxene (MgxFei_ySi206) type are rep- 
resent more than about 98% of silicate dust grainj^ in the ISM ( Kemper et al. 2005 : Min et al 



2007al ). where the protoplanetary dust grains are thought to originate. Olivine- type amorphous 
silicates show a broad triangular-shaped feature in the 10 /um region which peaks at 9.8 fim. Py- 
roxenes show a similar band to olivine in the 10 /im region, but its peak position is located at 
somewhat shorter wavelengths (~9.2^m). 

The broad features of amorphous silicates are less sensitive to the applied scattering the- 
ory (grain shape effects) than crystalline bands. It is thus not surprising that, apart from the 
size of the grains, not much information is available on the properties (e.g., shape, Mg-cqntent ) 
of the amorphous silic ate grains. For instance, most of the studies (e.g., Ivan Boekel et al.l (j2005l ). 



Bouwman et al 



(120081)') used the opti cal constants of iron-magnesium silicates with Fe/ (Mg+Fe)=0.5 



published by iDorschner et al.l (jl995l ) with Mie theory, assuming compact spheres for the grain 
shape. The aforementioned iron content of the silicates was used on the basis of cosmic element 
abundance constraints, and their higher mid-infrared opacities, compared to iron-free silicate grains. 
It is, however, s urprising that although protop lanetary dust grains are always regarded as porous 
aggregates (e.g., iHenning &: Stognienkd (jl99a )) a compact sphere model can fit the observed fea- 
tures rel atively well. From the analysis of the 10 ^m silicate absorption profile toward the Galactic 
Center (jMin et al.l l2007al ) concluded that the best fit can be obtained by using porous iron-free 
silicates. We used both iron-magnesium silicates with Fe/(Mg+Fe)=0.5 and iron-free silicates with 



^Here, we neglected all "featureless" dust species (e.g., iron and carbon) for which only weak constraints can be 
drawn from mid-infrared spectroscopy 
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Fe/(Mg+Fe)=0 and systematically tested the Mg-content of the amorphous grains and the scat- 
tering theory (i.e. grain shape effect). 

In contrast to the broad features of amorphous silicates, crystalline silicates show sharp and 
narrow features in the mid-infra red, which can be frequently seen in the spectra of both young 
and evolved stars (jHenningl l2009l ) . The analysis of the positions and the relative strength of these 
sharp features revealed that the radiating material s hould be a mixtu r e of fo r sterite (Mg2S i 04) an d 



enstatite (MgSiOs, see e.g., iMalfait et al.1 (|l998bl '): iBouwman et al.l (j200lh : iMeeus et al.l (1200l|)) 



These minerals are the magnesium-end members of the olivine and pyroxene solution series. Al- 
though crystalline silicates are usually minor dust components in protoplanetary disks compared 
to amorphous silicates, their sharp features can be seen in the spectrum in almost all cases. 

Studies of interplanetary dust particles show that these grains frequently contain large in- 
clusions of silica. Laboratory annealing experiments of amo rphous silicat e s also show that dur- 
ing the formation of forsterite, silica can be produced (e.g., iFabian et al.l ([20001)). Indeed silica 
has been found in the spe ctrum of youii g stars both in amorphous and in crystalline form (e.g., 
van Boekel et al.l (|2005l ) or lSargent et al.l (|2009al )). Silica shows a narrow, strong distinct features 
at ~9 iim and a broad, but also strong band at ~21 fim. The dust species together with the refer- 
ences of the applied optical constants are summarized in Table[2j Apart from the above-mentioned 
five dust species, we did not find any evidence for other abundant dust species in the Spitzer data. 

Three scattering theories were considered to calculate (MACs) from the optical constants, Mie 
scattering, continuous distribution of ellipsoids (CDE) and distribution of hollow spheres (DHS). 
These scattering theories are the most widely used methods to model mid-infrared spectra of young 
stars. We have two requirements for the computation method we wish to apply during the analysis. 
(1) The shapes and positions of the dust features in the Spitzer spectra should be reproduced as well 
as possible (2) The applied theory should also be valid outside of the Rayleigh limit. The reason for 
this second requirement is that we wanted to study the sizes of dust grains. In the strong bands at 
10 /xm one can already be outside of the Rayleigh limit for a micron-sized particle. The comparison 
of band position found in the spectra and those in calculated MACs rules out the Mie theory 
immediately (see, Table[3|. In Figures[I] and [2] we present the calculated absorption efficiencies of 
forsterite using differ ent scattering theories and compare them to laboratory measurements from 
Tamanai et al.l (j2009l ) . It can be seen, that dust band positions can be about as well matched with 
DHS as with CDE, since the calculat ed MACs do no t differ so much from each other than they do 
compared to Mie scattering (see also lMin et al.ll2003l ). Our second requirement, however, excludes 
CDE since it is strictly valid within the Rayleigh limit only. In the case of DHS both of our required 
conditions are fulfilled, and furthermore it is a fast computational method. We used, therefore, the 
DHS theory to calculate the MACs from the optical constants for our analysis. MACs of each dust 
species were calculated for three discrete grain sizes (0.1 ^m, 2.0 /im and 5.0 ^m). For forsterite 
and enstatite we used only two grain sizes (0.1 fim and 2.0 ^m) as we did not find any evidence for 
large (> 2.0 /xm) crystals. Silicate grains larger than ~5 /im are not considered, since they do not 
show feature in the studied wavelength range (5-35 fim). 
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In DHS one computes the scattering/absorption cross section of hollow spheres with a volume 
fraction / = Vtot/yvac, where Vtot is the total volume of the grain and Vvac is the volume of the vac- 
uum inclusion. The final MACs will then be an average over a whole distribution of hollow spheres 
with different values of /. It has already be shown that for crystalline silicates one should average 
over all possible values of / (fr om to 1. ) to g et the best agreement with the observed positions 
of crystalline bands (see e.g.,, iMin et al.l (120031 )). In Fig[3] we show the absorption efficiencies of 
amorphous silicates with olivine and pyroxene stoichiometry calculated using DHS theory. It can 
be seen that the higher the upper boundary for the hollow sphere distribution (/max) is chosen 
the broader the feature becomes. By increasing the value of /maxj the peak position of the feature 
shifts toward longer wavelengths. For the amorphous silicates we found that the best agreement 
with the observed spectra is obtained if one uses /max = 0.7. For the details, see Sec l4.2[ 



3.2. PAH band profiles 

All sources discussed in this paper show emission from PAHs. PAHs are also included in 
the spectral decomposition procedure in order to avoid systematic biases in the estimated dust 
parameters due to the PAH emission. PAH emission at 11.3 //m, 8.6 ^um and 12.7 /im can cause 
confusion in the estimated forsterite and silica content, respectively. In order to get the most 
realistic intensity profile for the observed PAH features, band profiles have been extracted from 
the spectra of sources with PAH emission only. These sources were HD34282, RR Tau, HD97048, 
HD135344B, HD141569 and HD169142. Five band profiles have been derived from the spectra of 
each source separately. We denote a set of profiles belonging to one source XI... X6, corresponding to 
HD34282....HD169142, respectively and we call the individual profiles after the central wavelength 
position as 6.2 ^um 7.7//m, 8.6 /im, 11.3 /im and 12.7 /im profiles (see FigureH]). The Xl-6.2 /im 
profile is therefore derived from HD34282 and its central wavelength is about 6.2 /xm. For further 
details of the derivation of the band profiles we refer to Acke et al. (in prep). 



3.3. Spectral analysis 

In order to analyze the dust composition in the disk atmosphere, the radiation of which domi- 
na tes the IRS spectrum , we used the two-layer temperature distribution (TLTD) method described 



m 



Juhasz et al.l (|2009l ). This method uses a multi-component continuum (star, inner rim, disk 
midplane) and it assumes that the region where the observed radiation originates (both optically 
thin and thick) has a distribution of temperatures instead of a single one. In this fitting method 
the observed flux density at a given frequency is given by 



N M rT ■ n 

I J^a.min 2 — qa 



i=l j = l •^Ta,r— 
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NP 
i=l 

where, N and M are the number of dust species and grain sizes, respectively. Np denotes the 
number of different PAH templates included in the fit. Kjj- is the mass absorption coefficient of 
the dust species i and grain size j. By{T) is the Planck function, qa is the power exponent of the 
temperature distribution and d is the distance to the source. The subscript a in the integration 
boundaries refers to the disk atmosphere. The continuum emission {Fy^cont) is given by 



F,,cont = ^B,{n) + Dl 



+ D2 



r ,max 

T 

-L m.min 



—B,{T)T—dT 

271" 2 — qm 

-^B4T)T—dT. 



(2) 



The first term on the right hand side describes the emission of the star, while the second and 
third terms describe the radiation of the inner rim and the disk midplane, respectively. The stellar 
emission, used for the fits, was not fitted during the mid-infrared spectral analysis, but it was 
derived from a separate fit to the UV-optical photometry from the literature. 

The assumptions (e.g., one single dust comp osition) used i n the TLTD method are not valid 
for an arbitrarily broad wavelength interval (see lJuhasz et al.l (|2009l )). Fitting the Spitzer IRS 
spectra to the total available wavelength interval (5.5-35 ^m) is already not reasonable. Therefore, 
we divided the Spitzer IRS wavelength range into two regions, 5.5-17 and 17-37 /im. These 
two wavelength intervals were fitted separately, although for the longer wavelengths we used the 
star and the rim emission which were fitted to the 5.5-17 ^m range. PAH templates were included 
in the fit only for the shorter wavelength interval. The final model in the 5.5-17 region was 
obtained using seven fits of each spectra. In the first fit only the XI set of PAH band profiles was 
used, in the second fit we used only the X2 profiles, etc. After the spectra were fitted with all six 
sets of PAH band profiles separately, we calculated the of the fit for the wavelength interval of 
the individual PAH bands. In the seventh fit we used a combination of PAH profiles taking the 
best-fit profile (with the lowest x^) for each band (e.g., Xl-6.2 /im, X5-7.7/xm, X6-8.6 /im, etc.). 
The final model for a given spectrum was chosen to be the one which gives the lowest global for 
the whole 5-17 //m interval. For the fits in the longer wavelength interval the rim contribution was 
not fitted, only the optically thin emission and the midplane component. We used the parameters 
for the rim which were derived from the fitting of the 10 /im region. 

To estimate the u ncertainties on the derived dust paramet ers we used a Monte Carlo type of 
error estimation (e.g.,, Ivan Boekel et al.ll2005l : iMin et al.ll2007al ). In this kind of error estimation, 
a normally distributed noise is added to the spectrum, scaling the width of the distribution to the 
simulated observational uncertainty in the flux value. Then the resulting spectrum is fitted. This 
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procedure was repeated 100 times. Then the standard deviation of the resulting mass fractions 
from the 100 fits wih be the uncertainty of the derived dust compositions. 

4. Results 

4.1. General summary of the fits 

The fitted dust composition for each source is presented in TableHlfTSl while the fits themselves 
are shown in Figure[5]|12[ The agreement between the observed spectra and our models are very 
good in general, with only a few exceptions. In three cases (HD35187, HD38120 and HD139614) 
a significant part of the in t^is long wavelength fits originates in the region between 30 fim and 
35 //m which is related to the problem with the 12th order of the long-high module (see Section [2?2]) . 
In the case of IID36917 our model has difficulties to match the observed spectrum longward of 
14 iim, which could be caused by the presence of a 16-19 //m PAH-band complex, which we did 
not take into account during the fitting. We believe, however, that these problems did not affect 
the main results of this paper. The reason is that the crystalline emission features in the long 
wavelength interval, which are investigated in details later on, are either very weak or completely 
missing in the spectra of these sources. 

For the rest of the sample differences between model and observation are usually at the per- 
centage level shortward of 17 //m and 5%-8% longward of 17 //m. The reduced values are, 
however, usually several tens in contrast to the expected value of about one for a good fit. We 
should keep in mind that the spectra analyzed in this paper have extremely high S/N (typically 
several hundreds). There are several effects which are negligible for lower S/N ratio spectra but 
become important for such extremely high S/N. The most important group of these effects is that 
during the calculation of the we took only the uncertainties of the Spitzer IRS spectra into 
account and we neglected all uncertainty related to our dust model. 

For instance, it is known that grain shape is an important parameter if dust grains are in the 
Rayleigh domain, especially for crystalline silicates. Protoplanetary dust grains are thought to have 
irregular shape where the calculation of the MACs from the optical constants are not straightfor- 
ward. Differences between MACs calculated by different scattering theories in the Rayleigh domain 
are much larger than a few percent (see Fig llSp . which is a typical discrepancy level in our fits. Our 
neglected uncertainty on the grain shape is also supported by the fact that the quality of the fit 
usually gets worse for spectra with higher crystallinity (see FigurelH]). 

Another source of uncertainty is the chemical composition of our dust model. We used the 
laboratory measurements of certain materials that are analogous to, but not necessarily the same 
as that in the astronomical environment. Slight differences in the chemical composition of the 
material (e.g., iron-content, Ca, Al or other ion inclusions) can already change the band profile 
at a percentage level (see Figure[3]). Even if the composition of the material is the same, their 
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band shapes are not necessarily identic a.1 (see for example the m eas urements of Mg-ric h amorphous 
silicates with pyroxene stoichiometry bv lDorschner et alj (|l995l ) and lJager et al.l (|2003l )). Therefore, 
we certainly found the abundant types of dust species and minerals, but we cannot claim that we 
found the exact material composition. All these types of uncertainties are real and are present in 
our data/analysis, however they cannot easily be measured and incorporated into the calculations. 



4.2. Amorphous silicates 



As a first step, we collected the most recent measurements of optical constants for amorphous 
silicates and teste d the iron content of th e dust grains together with the applied scattering theory 
(i.e. grain shape). iDorschner et al.l (|l995l ) published optical constants of glassy silicates with various 
iron content. Their measurements cover Fe/(Mg+Fe) ratios between and 0.6 for the pyroxene and 
between 0.5 and .6 for the olivine family. The other set of optical constants was determined by 
Jager et al.l (j2003l ) on amorphous silicates produced by the sol~gel method. In these experiments 
only iron-free silicates were measured. Iron-rich amorphous silicates were not tested in our analysis, 
since there are no laboratory measurements of optical constants of iron rich amorphous silicates 
for both olivine and pyroxene stoichiometry available. IDorschner et al.l (119951 ) measured iron-rich 
silicates only with olivine stoichiometry. We defined three mixtures of amorphous silicates optical 
constants to be tested. 



• AM IXl iron-free silicates w ith optical constants from I Jager et al.l (120031 ) for the olivine and 
from IDorschner et al.l ()1995l ) for pyroxene stoichiometry. 



AMI X2 iron-magnesium silicates (Fe/(Mg-|-Fe)=0.5) with optical constants from IDorschner et al 
(119951 ) for both olivine and pyroxene stoichiometry. 



AMIX3 iron-free silicates with optical constants from iJager et al.l (120031 ) for both olivine 
and pyroxene stoichiometry. 



We calculated the MACs from the optical constants using DHS theory for a grid of fn 
from (identical to Mie theory) to 1.0. 



values, 



In order to study the amorphous silicates in detail we selected three sources (HD36112, 
HD144432 and HD152404) where the mid-infrared dust features show the highest possible con- 
tribution from small amorphous silicate grains over any other optically thin emission. In other 
words, (i) emission from crystalline silicates should be the lowest possible, (2) emission of amor- 
phous grains should be dominated by small grains (< 1 //m), (3) contribution of PAH emission 
should be the lowest possible. We use the empirical "feature strength vs. shape" diagram of 
the W fim silicate feature for the selection (see Figure fTSl) . A third-order polynomial continuum 
is fitted to the 10 /im region fo r each spectrum and t he feature strength is then calculated as 
Fmax=l+(F°^'-F™'^*)/<F™'^* > (|van Boekel et ahlboosl ). and the feature shape is the ratio of the 
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continuum subtracted spectrum at 11.3 /im and 9.8 /im. Pristine W fim features lie in the bottom 
right corner of this diagram, while 10 /im complexes with the strongest contribution from large 
grains and crystalline silicates lie in the upper left corner. The selection criteria were Fmax > 3.2 
and Fii,3/F9.8 <0.71. 

Spectra of the selected sources were fitted using the MACs calculated from the AMIXl, AMIX2 
and AMIX3 mixtures separately for the amorphous silicate components. Fits were performed for 
each value of the fmax grid (i.e. different grain shapes). The MACs of the other dust components 
used in these fits were identical to those in the final fits and their mass fractions were allowed to 
vary. The results are summarized in Figure[T6l It can easily be seen that the AMIXl mixture of 
iron- free silicates gives always lower x^, than the other two mixtures. The increase in the if 
we used AMIX2 or AMIX3 instead of AMIXl was always far more than 3 a. We also noted that 
fits with AMIXl or AMIX3 (i.e. pure magnesium silicates) resulted in higher value of crystallinity 
compared to AMIX2. For the AMIXl mixture the best fit, i.e. lowest x^j is obtained if one used 
a DHS theory with fmax = 0.7, which means that the dust grains should be porous. It is also 
interesting to note, that using iron-magnesium silicates (AMIX2) less porous grains (i.e. lower 
values of fmax) usually give better fits than more porous grains. During the analysis of the whole 
sample of stars we used the iron- free amorphous silicates of AMIXl. 

The porous iron-free amorphous silicates can reproduce the observed amorphous silicates fea- 
ture better than compact iron-magnesium silicate grains. Although at first glance this suggests a 
possible lower iron-content of the amorphous silicates than it was previously thought, the conclu- 
sion is not straightforward. If dust grains contained the iron in the form of metal lic inclusions the 



prese nce of iron will not change the shape of the mid-infrared silicate bands (e.g., iQssenkopf et al 



19921 ). Moreover, amorphous silicates even with the same chemical compos ition can show slight 



differences in the band shape (see the discussion in e.g., I Jager et al.l l2003l ) . Although we used 
magnesium amorphous silicates for the spectral analysis we will not discuss the presence or lack of 
iron in further details. The fits show that the size of the amorphous grains varies between 0.4 fim 
and 5 //m indicating significant grain growth in the disk atmosphere. 



4.3. Crystalline silicates 



Sharp bands of crystalline silicates can be observed in almost all spectra indicating that thermal 
processing of protoplanetary dust has already taken place in these disks. The value of crystallinity 
ranges fr om ^1 % up to ^30 % , which is about the same range found in lower mass T Tauri stars 



(see e.g.. lBouwman et al.l (j2008l ) : iMeeus et al.l ((20091)). The size of the silicate crystals, derived from 
the short wavelength fits, is found to be significantly smaller than that of the amorphous grains. 
In the 17-35 /im wavelength interval it is hard to determine the size of the amorphous grains due 
to the very smooth mass absorption coefficient curve of the amorphous silicates. The size of the 
enstatite grains is usually larger than that of the forsterite grains. 
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As pointed out by lJuhasz et al.l (|2009l ). simple spectral analysis methods (like the TLTD 
method applied in this paper) are not capable of deriving dust parameter gradients in the disk 
directly. Some spatial information on the dust composition can, however, still be derived by com- 
paring the results of fits using the shorter (7-17 /im) and longer (17-35 /im) wavelength regions. 
By comparing the two regions we found that the ratio of the two crystalline dust species, enstatite 
and forsterite, changes. The shorter wavelength domain tends to show higher enstatite abundance 
compared to forsterite, than the longer wavelength bands. 



4.4. Correlations 

We searched for correlations between global system parameters and the properties of the dust 
grains in order to learn more about the mechanisms driving dust grain processing in general and to 
link dust and global system evolution. Such global parameters are the stellar parameters (luminos- 
ity, temperature), disk flaring, disk mass and the slope of the SED at sub-millimeter wavelengths. 
The disk flaring was empirically parameterized by the flux ratios at 24 fim and 8 /xm|. The disk 
mass and sub- millimeter slope were determined for only a part of the sample where sub- millimeter - 
millimeter observations were available. To investigate the strength of the correlation quantitatively, 
we calculated the Pearson's correlation coefficient. Pearson's correlation coefficient (r) measures 
the linear dependence between two variables (x, y) and it is calculated as 



where n is the number of data points. In order to investigate the significance of the Pearson's 
correlation coefficient, one can calculate the probability that the observed relation can be produced 
by a random distribution with the same sample size (n) 

^(^'") = ^rS^ /V - u^)^''-'^'^du (4) 
V^rP ( — j J\r\ 



(see, e.g... lTayloiill997l ) 



We did not find any correlation between stellar parameters (e.g.,, luminosity, temperature) and 
the dust parameters we derived. We want to note, however, that there is no large spread in these 
parameters within our sample, thus this conclusion should be taken with care. We also looked 
for possible correlation between the derived dust parameters from Spitzer IRS spectra and disk 
mass or slope of the SED at sub-millimeter wavelengths, for those sources where sub-millimeter 
measurements were available. No correlation was detected between fitted dust composition and 
disk mass or slope of the SED at sub-millimeter wavelengths. The lack of correlation suggests that 



^ We note that the exact choice of the wavelengths at which the flux ratio was taken does not affect our results. 
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the dust populations traced by mid-infrared and sub-millimeter wavelengths are not related to each 
other. 

A strong correlation was found between the mass-averaged grain size of the amorphous silicates 
and the flaring of the disk, which was empirically measured by the flux ratios at 24 ^um and 8 /xm. As 
can be seen in Figure fTHl that sources with flatter disks have larger grains in their disk atmosphere. 
We also investigated this correlati on for Group I and Gro up II sources. For the classification of 



the sources we used the criterion of Ivan Boekel et al.l (|2005l ) based on the near- infrared to infrared 



luminosity ratios and the IRAS 12-60 //m co l or (se e Fig llTp . This gives the same results as the 



original classification scheme bv lMeeus et al.l (|200ll ). but is somewhat easier to compute. Within 



Group II the Pearson's correlation coefficient is -0.65 with a probability that it is produced by a 
random distribution of 2 x 10~^. Although there are Group I sources which fit to the trend, defined 
by Group II sources, there are also several outliers. Investigation of the SEDs of these outlier 
Group I sources revealed that their disk structure is clearly different (due to e.g., the presence of 
a large inner hole or a gap) from that of those sources which fit to the trend. In these extreme 
Group I sources the calculated flux ratio is likely influenced by other disk parameters (e.g., size 
and location of the gap or inner hole) as well, and the flux ratio is, therefore, not a unique and 
independent tracer of disk flaring. Hence, we re-calculated the Pearson's correlation coefficient 
for the sub-sample of those Group I and Group II sources, which have 'usual' disk structure (i.e. 
without a gap or larger inner hole visible in the SED). The resulting correlation coefficient is -0.7 
with a probability of 6 x 10^^ that this correlation is produced by a random distribution. 

We did not detect any significant correlation between the value of crystallinity and any other 
global system parameter. However, we found that the enstatite over forsterite ratio, derived from 
the 5-17 /xm fits differ significantly in Group I and Group II sources, with higher values for Group 
II sources. The Kolmogorov-Smirnov test gives 2% probability that the enstatite-forsterite ratios 
of the two groups originate from the same distribution. The size of the crystals does not correlate 
with that of the amorphous grains, as one would expect if the crystals were produced by annealing. 
In contrast to the amorphous dust population, the grain size of the crystalline silicates does not 
correlate with the disk flaring. No evidence was found for correlation between amorphous grain 
size and crystallinity. 



5. Discussion 



5.1. Grain growth 



Initially the bulk of the dust in protoplanetary disks is present in sub-micron sized grains. 
Later on dust particles are thought to form larger aggrega tes and settle to the disk midplane (e.g., 
Dullemond &: Dominikll2004bl: ISchrapler &: Hennina 120041') . B v an alyzing the Spi t zer IR S spectra 



Bouwman et al 



(|200a) 



Meeus et al. 



(|2009|) and ISareent et al.l (|2009bl ^ found a 



of T Tauri stars, 

correlation between the average grain size of amorphous silicates and the disk flaring, parameterized 
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by the flux ratio between 30 ^um and 13 ^m. The question arises whether grain growth was in fact 
require d to produce this corre l ation or if dust setthng without coagulation can give the same 
result. iDullemond DominikI (|2008l ) studied this question, using two-dimensional axisymmetric 
disk models. Their results suggest that the effect of dust sedimentation alone (without coagulation) 
results in stronger, sharper 10 fim silicate features in flatter disks compared to the flared ones. As 
larger dust grains settle toward the disk midplane the disk becomes geometrically flatter, which 
can be observed as a decrease of the mid- to far-infrared flux in the SED compared to near-infrared 
fluxes. Since only the smallest grain population was left behind in the disk atmosphere, the emerging 
10 //m silicate feature should be strong and pointy similar to that observed in the ISM. We find the 
opposite trend, as the flux ratio between 24 /im and 8 fim decreases (flatter disks) t he average grain 



size of the amorphous si li cates increases. This is in agreement with the findings of lBouwman et al 



(|2008l ) and lMeeus et al.l (j2009l ) who also found the same trend in a sample of T Tauri stars. The 
observed correlation is, therefore, a proof that grain growth should have happened in these disks 
and that coagulation of dust particles and dust settling caused the fiattening of these disks. 

The observed correlation between disk flaring and grain growth can also be enhanced by a 
grain size gradient in the disk. Due to the decrease of the flaring index, the radial temperature 
(and therefore brightness) distribution in protoplanetary disks becomes steeper. This means that 
in flatter disks the mid-infrared emitting region is more concentrated toward the central star than 
in flared disks (e.g., iKessler-Silacci et al.l l2007). If the size of the grains increases toward the central 
star, flatter disks show even more evolved silicate features than they would in the lack of the grain 
size gradient. 

The found corr elation also supports the evolu tionary trend between Group I and Group II 
sources suggested by IDullemond &: DominikI (j2004al ) . While Group I sources have flared disk which 
is dominated by small grains, the average grain size in the disk of Group II sources is larger and 
dust settling causes the flattening of the disk. According to this prediction Group I sources should 
occupy the bottom right regions in Figure[18l while Group II sources should be located in the top 
left corner, which is indeed visible. We note again, that the outlier Group I sources in Figure[T8| 
which do not flt to the main trend, have different disk structure than the rest of the sample. 
Therefore, it is not straightforward to compare them to the main trend drawn by Group II sources. 

Although a tight correlation has been found between the disk flaring and the size of the 
amorphous grain population, no correlation was found between size of the crystals and grain size of 
amorphous silicates or disk flaring. A possible explanation for this can be if small silicate crystals 
are parts of larger aggregates. If the aggregate is very fluffy and the mass fraction of the crystalline 
constituents is low, crystalline bands will appear i n the mass abso rption coefficients as if crystals 
were isolated from the large amorphous aggregate (jMin et al.l 120081 ). 
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5.2. Silicate crystals 

It is an interesting and highly debated question how siUcate crystals form in protoplanetary 
disks. So far three scenarios have been pr oposed for the location of crystal formation in proto- 
planetary disks and the source of heating. iGaill (j2004l ) proposed that silicate crystals formed via 
thermal annealing and conde nsation in the ve r y imi er regions of protoplanetary disks due to the 
heating of viscous accretion. iHarker &: DeschI ([2002) proposed that crystalline silicates can form 
via thermal ann ealing at severa l AU f rom the central star if dust grains are heated by shocks. In 
a recent paper, 



Abraham et al 



(120091 ) reported that episodic crystal formation via annealing in 
the surface layers of protoplanetary disks during accretion outbursts can also be a possible way of 
producing crystalline silicates. 

From this study the following constraints, which should be explained by the model of crystal 
formation, can be identified; 



(1) Bands of crystalline silicates are seen in the 20-30 //m wavelength interval, suggesting that 
crystals are also located in low temperature (~100K) regions. 

(2) The forsterite-to-enstatite mass ratio changes with radius; its value is lower in the shorter 
wavelength fits (inner disk) compared to that in the longer wavelengths (outer disk, see 
FigurettH) . 



(3) Crystallinity does not correlate with any global parameter of the system (e.g., stellar 
parameters, disk mass, disk flaring). 

(4) The size of the crystals is smaller than that of the amorphous grains. 

(5) The grain size of crystalline silicates does not correlate with that of the amorphous grains. 

(6) No strong evidence has been found for the correlation of crystallinity and amorphous grain 
size. 

(7) The size of the enstatite grains is usually larger than that of the forsterite crystals. 



All three n i odels of crystal formation are consistent with (1). In the model of iGaill (j2004l ) and 
Abraham et al. ( 2009 ) crystals are formed in the inner disk (<1-2AU) and then transported to 
the outer, colder regions by large-scale r nixing . In th e case of shock heating crystals are produced 
in the outer regions. In the models of iGaill (j2004l ) crystal formation occurs under equilibrium 
conditions. Under chemical equilibrium the crystal population is dominated by enstatite, assuming 
solar element abundances. Forsterite can only be the dominant cr ystalline com ponent in the very 
inner disk due to its somewhat higher stability limit than enstatite. iGaill (|2004l ). therefore, predicts 
decreasing abundance of forsterite relative to enstatite as a function of radius in protoplanetary 
disks. Our finding (2) contradicts this prediction, suggesting that the crystals we see in the spectra 
of HAe stars were not formed under equilibrium conditions. If crystals were produced by shock 
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heating the time is probably too short t o achieve chemical e quihbrium. This has already been proved 
by laboratory annealing experiments ([Fabian et alj 120001 ) . These experiments demonstrated that 
annealing of small (0.1-1 ^m sized) porous silicate grains with olivine and pyroxene stoichiometry 
results in the formation of forsterite regardless of the starting stoichiometry. These experiments 
also showed that the formation of forsterite, in the case of these small particles, is about an order 
of magnitude faster if the starting stoichiometry of the amorphous silicates is olivine instead of 
pyroxene. In the case of shock heating, where the annealing time is very short (less than an hour), 
one would expect the transformation of amorphous silicates with olivine stoichiometry to forsterite. 
Thus, the resulting crystal population will be dominated by forsterite while in the a morphous phase 



the ab undance of grains with pyroxene stoichiometry should increase. The model of lAbraham et al 



([20091) does not have a clear prediction for the forsterite-to-enstatite ratio. Depending on the 
duration of the outburst and on the time a dust grain spent in t he di sk atmosphere the mass 
of the produced enstatite crystals can vary. [Bouwman et al.[ ([2008[ ) and [Meeus et al.l ([2009 ) also 
speculated on the formation of enstatite. Our findings support their conclusion that the formation 
of enstatite probably occurs at the very inner regions of protoplanetary disks. The temperature can 
be high enough to reach chemical equilibrium during the crystal formation, favoring the formation 
of enstatite instead of forsterite. Moreover, the higher densities of the inner disk compared to the 
outer regions can result in frequent collisions betw een dust particles mak ing the dust grains more 
compact which also favors the enstatite formation ([Bouwman et al.ll2008l ) 



In the model of lGaill ([2004 ) the required heating power for the crystallization is provided by the 
viscous heating of accretion in the disk midplane. This model, thus, predicts a correlation between 
the mass of the crystals and the accretion rate. It is, however, not obvious if this scenario predicts 
any correlation between the observed crystallinity and any of the global parameters of the system. In 
this scenario, crystals form in the disk midplane, thus the amount of crystals, that can be observed 
in the disk atmosphere, depends on the efficiency of both the radial and the vertical mixing in the 
disk, which are basically unknown factors. Since shocks in protoplanetary disks can be produ ced by 
several different mechanisms (e.g., bow shock of a planetary embryo, see iDesch et al.l ([20051 )) and 
they can act at the same time, this model does not predict a strong correlation between crystallinity 
and global parameters of the system. The episodic crystal formation scenario, in agreement with 
(3), predicts even less correlation between crystallinity and any parameter of the star or the disk. 
Since outbursts, which provide the heating power for the crystallization, happen randomly and the 
total number of outburst a system goes through is also unpredictable, this scenario predicts a large 
scatter in the value of crystallinity even within stars with similar parameters. 

Our findings (4) and (5) can be consis tent with b oth shock heating and midplane crystallization. 
If crystals were formed via condensation ([Gailll2004l ) one would not expect any correlation between 
the sizes of the crystals and the sizes of the initially amorphous material. This is due to the fact 
that the evaporation and re-condensation of dust grains erases all correlation between the size 
crystals and the parent amorphous grains. If crystal formation occurs via annealing one would 
expect, however, a correlation between the sizes of the starting and the resulting materials. This 
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correlation can only be erased if large amorphous aggregates are disrupted before or during the 
annealing process. During interactions with a shock, dust grains are not only annealed but larger 
aggregates are also likely to be disrupted. The details of the aggregate disruption depend on the 
physical conditions in the shock, which can vary significantly in protoplanetary disks. Thus, shock 
heating predicts weak or no correlation between the size of the crystals and the amorphous grains. 
In episodic crystal formation, however, it is not obvious how amorphous grains can be annealed 
and at the same time the grain size can also be changed. While shock heating and midplane 
crystallization naturally explain the lack of correlation between crystallinity and amorphous grain 
size (6), in case of episodic surface crystallization still remains to be explained. 



5.3. Quality of the dust model 



One of the main questions, we tried to answer in this study, is whether the dust model in 
protoplanetary disks we derived from lower S/N data is sufficient to fit the spectra with high S/N 
ratios. In other words, will spectra with higher S/N reveal new dust components or does the 
current dust model, previously used to analyze spectra of protoplanetary disks, already contain all 
the main dust species which show features in the mid-infrared domain? In our analysis we did not 
find any concluding evidence for new dust species in our high S/N spectra. However, we found 
small differences between the peak positions seen in the Spitzer IRS spectra and in our calculated 
absorption coefficients. 

The spectrum of everal sources (e.g.,, HD190073, IID35929, IID244604) show a smooth 10 ^um 
silicate feature, which does not show any substructure apart from a strong 9.3-9.4 //m peak. It is 
also interesting to note that in the spectra of these sources one can find several small peaks in the 
13-16 /xm wavelength interval, suggesting the presence of pyroxene-type crystals in the spectrum. 
Ortho- or clino-enstatite cannot reproduce the strong 9.3 ^m peak and the smooth substructure- 
less 10 /im feature at the same time. We derived the crystalline pyroxene MACs from the observed 
spectra. This has been done by subtracting all other fitted components (continuum, amorphous 
silicate, forsterite, silica) from the observed spectrum and dividing the resulting curve by the 
underlying source function. We derived the MAC curve from the spectra of four sources (HD190073, 
IID244604, HD35929, IID179218), where the 9.3^m feature was the strongest. The average curve 
is presented in Figure[20l We compared the der ived MAC curve to that of ortho- and clino- 
enstatite and an iron bearing crystalline pyroxene (jChihara et al.ll2002l ). Figure[20] shows that the 
de rived MAC curve is similar to that of pyroxene crystals co ntaining 10 % iron ( denoted by En90 



m 



Chihara et al.l 12002) . A similar comparison was done by iBowey et al.l (j2007l ) , who compared 



the peak positions of pyroxene crystals with diffe r ent ir on content to the mid-infrared spectra of 
IID104237 and IID179218 from Ivan Boekel et al.l jlOod ). On the basis of this comparison these 
authors suggested that the pyroxene crystals around these stars lik e ly cont ain 10 %-25 % iron. 
By analyzing the Spitzer IRS spectra of T Tauri stars I Sargent et al.l (|2009al ) also arrived to the 
conclusion that crystalline pyroxenes are unlikely to be pure magnesiasilicates. This implies that 
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crystalline silicates in protoplanetary disks are not necessary iron- free as it is usually assumed. 

Almost all sources show a 16 /im feature usually associated with crystalline forsterite. Although 
our dust model can fit the other forsterite bands in the 7-17 ^m interval quite well, the 16 //m 
band was never fitted with the same quality as the rest of the spectrum in this interval. In 
order to investigate this phenomenon in detail, we derived the forsterite absorption coefficients 
from the spectrum of HD100546, which shows the strongest contribution from forsterite in the 7- 
17 //m wavelength range. The derivation process was similar to that applied to derive the enstatite 
absorption coefficients. We subtracted all other contributions (continuum, PAHs, silica, amorphous 
silicates, enstatite) from the observed spectrum. The resulting curve is presented in Figure[2TJ We 
also compared the derived MACs to the calculated ones which have been used for the fitting and 



those measured in the laboratory. For the comparison we used the measurements of lTamanai et al 



(|2009l ) where the absorption coefficients were measured on free-flying particles instead of embedded 
in a KBr pellet. It can be seen that the position of the 16 fim band is more consistent with the 
MACs measured on free-flying particles. On the other hand, the peak position of the 10.0 fim and 
the 11.3 fim bands are reproduced better in the MACs calculated using DHS theory. In general the 
calculated MACs matched the observed bands better in the 17-35 /im wavelength interval than in 
the 5-17 /im region. However, in two cases (HD244604 and HD203024) the forsterite bands also 
longward of 22 //m seem to be shifted toward shorter wavelengths, similarly to the 16 /im feature. 
The reason for the difference between the peak position of the calculated and the observed forsterite 
bands is probably twofold. Intrinsic differences in the exact shape of the forsterite crystals result 
in shifts of the band position of about 0.1 um or eveii more , depending on the band and the shape 



of the particles (iTamanai et al.ll2009l : iMutschke et al.ll2009l ). Another possible reason can be found 



in the scattering theory w e used to calculate the MACs from the optical constants. As pointed out 
by Mutschke et al. (j2009l ) calculating the absorption cross section along each crystallographic axis 



independently and then taking the average can result in the wrong band positions for anisotropic 
irregular particles. 



6. Conclusions 

In this paper we analyzed the Spitzer IRS spectra of a comprehensive sample of Herbig Ae 
stars. The spectra have very high S/N ratios (usually of the order of several hundred), which 
allowed us not only to study dust evolution in the disks around HAe stars, but also to investigate 
if the high S/N spectra reveal any new dust component. From this study we concluded, that: 

• On the basis of the analyzed spectra, the major constituents of protoplanetary dust around 
Herbig Ae/Be stars are identifled as amorphous silicates with olivine and pyroxene stoichiom- 
etry, crystalline forsterite and enstatite and silica. No strong evidence for new dust species has 
been found. However, we found slight deviations in the peak positions of the crystalline band 
seen in the Spitzer IRS spectra and in our calculated MACs. The position of the 16 //m band 
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of crystalline forsterite always peaks at shorter wavelengths than m our calculated MACs. 
In some cases the 24 /um band of forsterite is also shifted toward shorter wavelengths in the 
spectra compared to the MACs. The 8-14 fim region of several spectra shows only one narrow 
peak at ~9.4 /im on top of the broad, smooth amorphous feature. We found that the 9.4 fim 
peak can most likely be associated with crystalline pyroxene with 10% iron content. 

A tight correlation was found between the average grain size of the amorphous silicates and 
the flux ratios between 24 ^m and 8 /x m, used as a proxy o f dis k flaring. Th i s find ing is 



consistent with what has been found bv lBouwman et al.l (|2008l ) and lMeeus et al.l (j2009l ) for T 



Tauri stars. The importance of this correlation is twofold. First, it is a strong observational 
proof that grain growth is required to produce the observed diversity in grain sizes derived 
from mid-infrared spectroscopy. Second, it also a strict proof that coagulation of dust grains 
and the accompanying sedimentation are responsible for the flattening of the disks observed 
as the decrease of the mid- to far- infrared radiation in disks. 

Those Group I sources which do not fit to the main trend between grain size and disk flaring 
were found to have different disk structure than the rest of the sample. The outlier Group I 
sources have large inner holes, or gaps which are visible in the SED. In these sources the 
calculated flux ratio between 24 /im and 8 fxm does not only measure the flaring index, but 
also it is strongly affected by other disk parameters (e.g., size and location of the inner hole 
or gap). 

We compared the predictions of currently existing theories for crystallization to the derived 
dust parameters and the potential correlations of these parameters with each other and global 
system parameters. We found that none of the three investigated crystallization mechanisms 
can alone reproduce all the observed correlations at the same time. It seems, therefore, very 
likely that all three processes may play at least some role in the evolution of protoplanetary 
dust particles. 



Crystallization in the disk midplane by accretion heating (iGaill l2004l ) predicts increasing 
enstatite-to-forsterite ratios with radius, in contrast to what is observed in the spectra of 
HAe stars. This scenario also has problems reproducing the observed lack of correlation of 
crystallinity with any stellar or disk parameters. If the accretion heating was responsible for 
the heating of the amorphous grains to crystallize them, the produced mass of crystals should 
correlate with the accretion rate, which is correlated with the stellar mass. Therefore, we 
should have observed increasing crystallinity with increasing stellar mass, which is not the 
case. 

In crystal formation by shock heating in the outer disk (at a few tens of AUs) amorphous 
dust grains are heated above the annealing temperature only for a very short time (coupl e 



of minutes to hours) during which the crystallization should occur (jHarker &: DeschI |2002| ) . 
Over such a short time-scale chemical equilibrium cannot be achieved and the resulting crys- 
tal product will be forsterite independently of the starting stoichiometry of the amorphous 
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particles. This can explain the observed dominance of forsterite in the outer disk, but it 
cannot explain the existence of enstatite at the same location. The formation of enstatite via 
annealing from forsterite is a slow process and it requires too long time to be locally produced 
by shock heating. 

T he recently sugge s ted ep isodic crystal formation in the surface layers of protoplanetary disks 



bv [Abraham et al.l (|2009l ) can naturally explain the observed diversity in the dust parameters 
and the lack of correlation of crystallinity with basically all the global parameters of the 
system. However, the observed lack of correlation between crystallinity and the size of the 
amorphous grains and the smaller size of crystals compared to amorphous grains still remains 
to be explained in this framework. 



We thank the anonymous referee for the careful review of the manuscript that helped to 
improve the paper. 
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Table 1. Source list of the sample and their coordinates. 



Name 


Ra (J2000) 


Dec (J2000) 


Associated 


Spec, type 


Group 


Dist 




ro ' "1 

L J 


ro ' "1 

L 1 


Region 








AB Aur 


4 55 45.84" 


+30 33 04.5" 


L1519 


AOVe+sh'^ 


la 


144.0" 


HD31648 


4 58 46.26'^ 


+29 50 37.2" 




A3ep+sh'= 


Ila 


131.0" 


HD34282 


5 16 00.48" 


-9 48 35.4" 




AOe^ 


lb 


164.0" 


HD35187 


5 24 01. 17'^ 


+24 57 37.8" 


L1559 


A2/3IV/Ve'^ 


Ila 


150.0" 


HD35929 


5 27 42.79"* 


-8 19 38.4" 


Orion OBlc 


FOIIIe'= 


Ila 


510.0^ 


HD36112 


5 30 27.53" 


+25 19 57.3" 




A5IVe^ 


la 


204.0" 


HD244604 


5 31 57.24^ 


11 17 41.46'' 


Orion OBla 


AO*^ 


Ila 


336.0^ 


HD36917 


5 34 46.99'' 


-5 34 14.59'' 


Orion OBlc 


AOV'' 


lib 


510.0- 


HD37258 


5 36 59.25^ 


-6 9 16.42'' 


Orion OBlc 


A2V'^ 


Ila 


510.0'= 


BF Ori 


5 37 13.26" 


-6 35 00.6" 


Orion OBlc 


A5-6IIIe'= 


Ila 


510.0^' 


HD37357 


5 37 47.09^ 


-6 42 30.17'' 


Orion OBlc 


AOV^ 


Ila 


480.0' 


HD37411 


5 38 14.51^ 


-5 25 13.26'' 


Orion OBlc 


Bgv'' 


lib 


510.0' 


RR Tau 


5 39 30.52^ 


26 22 26.98'' 


L153 


AOIVev'^ 


lib 


800.0" 


HD37806 


5 41 02.29" 


-2 43 00.7" 


Orion OBlb 


A2Vpe^ 


Ila 


470.0*^ 


HD38120 


5 43 11.89° 


-4 59 49.9" 


Orion OBlc 


B9/ 


la 


422.0" 


HD250550 


6 01 59.99" 


+ 16 30 56.8" 


Gem OBI 


B4-5IIIe'= 


la 


700.0' 


HD259431 


6 33 05.19" 


+ 10 19 20.0" 


Mon OBI 


BlVe'^ 


la 


290.0" 


V590 Mon 


6 40 44.64^ 


+9 48 2.1'' 


NGC 2264 


B8ep+shs 


la 


800.0™ 


HD50138 


6 51 33.40" 


-6 57 59.5" 




B5V[e]^ 


Ila 


289.0" 


NX Pup 


7 19 28.26" 


-44 35 11.4" 




A9-F0Ve^ 


Ila 


450.0^= 


HD58647 


7 25 56.10" 


-14 10 43.5" 




B9IVe'= 


Ila 


277.0" 


HD72106 


8 29 34.90" 


-38 36 21.2" 




AO-'^ 


Ila 


288.0" 


HD85567 


9 50 28.55" 


-60 58 03.0" 




B7-8Ve^ 


Ila 


1000.0" 


HD95881 


11 1 57.64'' 


-71 30 48.4'' 


Sco OB2-4 


AlIII"^ 


Ila 


118.0^ 


HD97048 


11 08 03.37" 


-77 39 17.5" 


Ced 111 


B9-A0ep+sh^ 


lb 


175.0" 


HD98922 


11 22 31.68" 


-53 22 11.5" 




Bgye-^' 


Ila 


1042.0" 


HD100453 


11 33 05.61" 


-54 19 28.5" 




A9V'' 


lb 


111.0" 


HD100546 


11 33 25.51" 


-70 11 41.2" 




B9Vne'= 


la 


103.0" 


HD101412 


11 39 44.46'' 


-60 10 27.84'' 


Sco OB2-4 


AO-'^ 


Ila 


160.0^ 


HD104237 


12 00 05.20" 


-78 11 34.5" 


Cha III 


A4IVe+sh^ 


Ila 


116.0" 


SS73 44 


15 3 23.81'' 


-63 22 58.87'' 


G317-4 


Beqs 


la 


1000.0" 


HD135344 


15 15 48.94'' 


-37 8 55.86'' 


Sco OB2-3 


F4V'' 


lb 


140.0^' 


HD139614 


15 40 46.39'' 


-42 29 53.63'' 


Sco OB2-3 


A7Ve'* 


la 


140.0^' 


HD141569 


15 49 57.76" 


-03 55 16.2" 


L169 


B9.5Ve'= 


lib 


99.0" 


HD142666 


15 56 40.02'' 


-22 1 40.01'' 


Sco OB2-2 


A8V'^ 


Ila 


145.0^ 


HD142527 


15 56 41.90" 


-42 19 23.1" 


Sco OB2-3 


F7IIIe^ 


la 


198.0" 


HD144432 


16 06 57.96" 


-27 43 09.6" 


Sco OB2-2 


A5Ve'= 


Ila 


253.0" 



Table 1 — Continued 



Name 


Ra (J2000) 


Dec (J2000) 


Associated 


Spec, type 


Group 


Dist 






r 1 


r ' "] 


Region 






[pc] 


HD144668 


16 


08 34.29° 


-39 06 18. r 


Lupus III 


A7IVe'^ 


Ila 


208.0'^ 


HD145263 


16 


10 55. 11'^ 


-25 31 21.4" 


Sco OB2-2 


FOV'' 


Ila 


116.0" 


Wray 15-1484 


16 


27 15.11'' 


-48 39 26.83'' 




BQS 


la 


750.0° 


HD150193 


16 


40 17.93^^ 


-23 53 45.0^^ 


Sco OB2-2 


AlYe'= 


Ila 


150.0" 


HD152404 


16 


54 44.86'^ 


-36 53 18.3-^ 


Upper Sco 


F5-fF5IVe'= 


Ila 


145.0" 


KK Oph 


17 


10 8.110'' 


-27 15 19.01'' 


Sco OB2-2 


A8Vev'= 


Ila 


160.0'= 


51 Oph 


17 


31 24.95'^ 


-23 57 45.3-^ 




B9.5Vne'= 


Ila 


131.0" 


HD163296 


17 


56 21. 29'^ 


-21 57 21. 5'^ 




AlVe'^ 


Ila 


122.0" 


HD169142 


18 


24 29.79'' 


-29 46 49.22'' 


Sco 0B2-1 


B9Ve'^ 


lb 


145.0^ 


VV Ser 


18 


28 47.86'' 


+0 8 39.77'' 


Ser Rl 


AOVevp^ 


Ila 


330.0^ 


T CrA 


19 


1 58.78'' 


-36 57 49.89'' 


NGC 6729 


B9e9 


la 


130.0^= 


HD179218 


19 


11 11. 25'^ 


+15 47 15. 8-^ 


L693 


B9e'= 


la 


244.0" 


WW Vul 


19 


25 58.75'' 


21 12 31.35'' 


Vul Rl 


A2IVe'= 


Ila 


440.0^' 


HD190073 


20 


03 02. 5r 


+5 44 16.7'^ 




AOIVepsh'^ 


Ila 


290.0^' 


LkHa 224 


20 


20 29.36'' 


41 21 28.44'' 


Cyg Rl 


B2/3e9 


la 


980.0'= 


HD203024 


21 


16 2.990'' 


68 54 52.08'' 


Cep R2 


B8.5V' 


la 


450.0' 



References. 



(d) 
(hi 



ESA 



Malfait et all (ll998aD (e 



Mannings fc Barlow! ( 1998 ). (i) iKun et al.i (120001) . (\) 



( 20061). m 



Wade et al 



9971). (b'l 



Mora et al 



Cutri et al, 
4200 1"), (f 



mm. (c) 



Vieira et al 



llYie 



20071) . r m) iHernandez et all (|2004l ). (n) 



van den Ancker et al. 
( 20031). (g ll The et al 



Acke et al 



1998), 



1994), 



( 2005f). fk) iManoia et al 



Carkner et all (|l998l) (o) 



Acke fc van den Ancked (|2004l ). (p) iHillenbrandl (|l997[ ) 
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Table 2. Overview of dust species used. For each component we specify its amorphous vs. 
crystaUine structure, chemical composition, shape and reference to the laboratory measurements 

of the optical constants. 



# 


Species 


State 


Chemical 
formula 


Shape 


Ref 


1 


Amorphous silicate 
(Olivine stoichiometry) 


A 


Mg2Si04 


Hollow Sphere 


(1) 


2 


Amorphous silicate 
(Olivine stoichiometry) 


A 


MgFeSi04 


Hollow Sphere 


(2) 


3 


Amorphous silicate 
(Pyroxene stoichiometry) 


A 


MgSiOg 


Hollow Sphere 


(1,2) 


4 


Amorphous silicate 
(Pyroxene stoichiometry) 


A 


MgFeSiaOe 


Hollow Sphere 


(2) 


5 


Forsterite 


C 


Mg2Si04 


Hollow Sphere 


(3) 


6 


Ortho-Enstatite 


C 


MgSiOs 


Hollow Sphere 


(4) 


7 


Silica 


A 


Si02 


Hollow Sphere 


(5) 



References. — (1) 
Sogawa et aP |2006l ): (4) 



Jager et al. 



Jager et al. 



(l2003h : (2 ) borschner et all dlQQSh (3) 
(|l998h : (51 iHenning fc Mutschkel dlQQj 
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Table 3. Summary of the strongest, most frequently obs erved features i n the Spit zer spectra o f 



HAe st ars. The identifications were done on the basis of IChihara et al.l (120021 ) and lKoike et al 



(j2003l ) for enstatite and f orsterite, respe c tively . The abbreviations in the case of enstatite are 



identical to those used in IChihara et al.l (j2002l ) (CEn - clino-enstatite, OEn - ortho-enstatite, 



En90 - crystalline pyroxene with 10% iron content). Amin and Amax are the shortest and longest 
wavelengths at which the feature is observed. < Ac > is the average peak position. 



< Ac > 


Amin 


Amax 


Nr. of sources 


Identification 


6.26 


6.22 


6.28 


31 


PAH 


7.89 


7.81 


8.03 


21 


PAH 


8.22 


8.08 


8.45 


21 


Silica 


8.67 


8.57 


8.77 


24 


PAH 


9.34 


9.17 


9.49 


26 


Enstatite(CEn,OEn,En90) 


9.91 


9.69 


10.1 


22 


Forsterite, Enstatite(CEn,OEn,En90) 


11.24 


11.14 


11.43 


48 


Forsterite, PAH 


12.68 


12.60 


12.84 


29 


PAH, Sihca 


13.76 


13.54 


13.90 


20 


Enstatite(CEn,OEn,En90) 


14.43 


14.38 


14.46 


12 


Enstatite(CEn,OEn,En90) 


14.63 


14.57 


14.71 


18 


Enstatite(CEn,En90) 


15.56 


15.35 


15.72 


23 


Enstatite(CEn,OEn,En90) 


16.12 


16.01 


16.33 


29 


Forsterite 


23.74 


23.00 


24.07 


21 


Forsterite, Enstatite(OEn,En90) 


27.47 


27.00 


28.20 


22 


Forsterite, Enstatite(Cen,OEn,En90) 


33.79 


33.39 


34.13 


19 


Forsterite, Enstatite(Cen,Oen,En90) 
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Table 4. Best fit dust parameters and the reduced of the fit (5-17 //m). Listed are the derived 
mass fractions of each dust component in % and the fitted PAH fluxes in Jy. The tabulated mass 

fractions (PAH fluxes) and their uncertainties are rounded to the nearest hundredths 
(thousandths) pl ace. In the case of the spectra where we used TIMMI2 spectra from 



van Boekel et al.l (|2005l ) due to the lack of Spitzer short-high module we excluded the 6.2 /xm PAH 
band from the fitting, which is indicated by a dash (-). 



AB Aur HD31648 HD35187 HD35929 HD36112 HD244604 HD36917 











80.83 




110.38 




41.94 




19.26 


54.24 




82.45 




89.55 


Am. Ol. 





1 /jm 


48 


gg + 6.70 










Q^ + 0.11 

'J^-0.01 





QQ+0.00 


49.40li-79 





00+°°° 
'-"^-0.00 





QQ+0.00 
'^'-'-0.00 


Am. Ol. 


2 


/tm 


31 


^3-4.64 


41 


6011:^3 


71 


Q5+0-81 

'''-'-0.78 


16 


41+2-74 
^^-2.43 


13.90tij9 





QQ+0.00 

'J'^-o.oo 





00+°°° 
^*^-0.00 


Am. Ol. 


5 


/jm 





09+0.04 
^^^-0.02 


22 


00+3.50 

*'^-3.94 





OQ+0.16 

'J'-'-o.oo 


47 


00+4. 09 
0»-3.93 


02+° °'^ 
'-'■'J^-0.02 


64 


,0+0.44 

^■^-0.43 


86 


o,+0.47 
^^-0.51 


Am. Py. 





1 /^m 


13 


00+2.59 
■^"-2.49 





QQ+0.00 

^'-'-0.00 





00+0.00 
'^'-'-0.00 





00+0.00 
'-"^-0.00 


24 54+0-79 





QQ+0.00 

^'^-O.OO 





nq+0.15 
"''-0.03 


Am. Py. 


2 


lira 





'-'^-0.02 


16 


gg+0.62 
"°-0.S6 


5 


on+0.70 

*''J-0.92 





n'3+0.76 
^'-'-0.03 


9 fi5+0.67 
^■°°-0.67 


8 


-15+0.40 
^"-'-0.29 





QQ + 0.00 
UU_Q 00 


Am. Py. 


5 


lira 





QQ+0.00 

^'^-O.OO 





= 0+2. 01 
"^•^-0.50 


1 


90+1-3? 


14 


40+1.75 

^'^-1.75 


00+° °° 

^■'^^-o.oo 


17 


QQ+0.61 

^^^-0.74 





QQ + 0.00 

'^^^-O.OO 


Fors. 





1 lira 





/iq+0.30 
^'-'-0.28 


1 


94+0 06 

^^^-0.06 


1 


00+0.06 
^■^-0.06 





'"'-0.07 


9 ic:+0.03 
^■^'-'-0.03 


1 


QK + 0.02 

"^-0.02 





QQ+0.00 

'^^-o.oo 


Fors. 


2 


lira 





vo+0.67 
'*'-0.64 





07+0.24 
-0.24 





00+0.17 
^^-0.02 





Qg+0.28 
'^"-0.06 


00+° °° 
'-'-^'-'-0.00 





QQ+0.00 

"^'^-O.OO 


6 


40+0.27 

^■'-0.25 


Enst. 





1 lira 





fio+o.es 


2 


1C+0.09 
^^-0.12 





01 +0.06 
"-^-0.05 


2 


11+0.12 
-^-^-0.10 


1 Oq+0.07 

-^■^■-'-0.07 


1 


00+0. 03 
°O-0.04 





QQ + 0.00 

'J'^-o.oo 


Enst 2 


lira 





in+0.65 
^'J-o.io 


3 


"^-0.21 


3 


n^+0.17 
'^^-0.20 


9 


00+0.53 
■-'°-0.45 


2 59+°-20 

Z.O»_Q 22 


4 


cq+O.lO 
'-'■'-0.09 





QQ+0.00 

""^-O-OO 


Silica 


1 lira 





1 0+0.46 
^^-0.11 


1 


ni+o.09 

'^^-0.06 





49+0.02 
^^-0.02 





cfi+O.OS 
'^°-0.06 


in+°-°3 

'-'■^^'-'-0.02 


1 


tri+0.02 
■^^-0.02 





nr;+0.02 
'^'-'-0.03 


Silica 2 


lira 


1 


1 7+0.76 
-^'-0.76 


2 


74+0.15 
'^-0.17 





09+0.06 
'-'^-0.02 





00+0.00 

'J'J-o.oo 


49+0.12 
■^■^^-0.11 





Qq+0.04 
='■'-0.03 





QQ + 0.02 

'J'J-o.oo 


Silica 5 


lira 


3 


04+1.92 
■^^-1.70 


7 


09+0-63 
^^^-0.68 


15 


9n+0.43 

^^^-0.43 


9 


nt; + 0.70 
'-'^-0.62 


00+°-°° 
'J-'J^-o.oo 





Q4+0.11 
'J^-0.03 


6 


00+0.27 

t'*'-0.29 


PAH 6.2 


3.089ing 




^62+''-°°5 


000+°-'"'° 

U.UUU_Q QQQ 


901+0.003 

U.ZOi_Q 


077+0.002 


990+0.003 

U.ZZO_Q 003 


PAH 7.7 




'J-^^^-o.oii 


625+''-°°'' 

U.DZO_Q QQg 


001+°-°°! 

U.UU1_Q 


n 045+°-°°^ 


n n4s+°-°°2 
U.U4S_Q 002 


409+°-°°2 

U.4UM_Q 002 


PAH 8.6 


2.6391E;-525 




400+0.007 


000+°-°°° 

U.UUU_Q QQQ 


n Qotr+0.009 

U.UJO_g QQg 


-1^:5+0.002 

U.100_Q 002 


,07+0.002 

U.1»/_Q.002 


PAH 11.3 


q .90 + 0.342 


0.428t0;«l* 


040+0.014 


007 + 0.003 


000+° °°° 

U.UUU_Q QQQ 


0.0001°;°°° 


n i;a7+0 004 
U.&»/_Q 005 


PAH 12.7 


U.Dli_g JJ3 


0.098tO:«f, 


1 ,0 + 0.006 

>J.iJ-Z_g QQg 


ni 1+0.001 

U.Uii_Q 


046+° °°^^ 

U.U4D_Q QQg 


0.089t°;°°^ 


203+°-°°^ 

U.ZUO_Q 003 
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Table 5. Same as TableU 



HD37258 BF Ori HD37357 HD37806 HD38120 HD250550 V590 Mon 











45.63 




68.54 




31.36 


129.77 




119.98 






39.36 




194.36 


Am. Ol. 


0. 


1 fim 


0, 


00+°'°° 


12, 


ir,+ 1.95 

■-^^-2.28 


0, 


00+°°° 
■'^'-'-0.00 


00+° °° 


25.4ll^-|6 




4, 


5-1+0.43 
■''i-o.si 


0, 


00+°°° 
■'^'-'-0.00 


Am. Ol. 


2, 


fim 


47, 


09+0.54 
■"^-0.58 


38, 




61, 


70+0.50 
"^-0.52 


9 91+19.27 
^•^1-2.14 


59.16tt:" 


33, 


97+1-44 
-1.23 


27, 


79+0.69 
■ '^-0.75 


Am. Ol. 


5, 


,0 fim 


0, 


Qn+0.12 
■'^'^-0.00 


0, 


QQ+O.OO 


0, 


00+°°° 

■'^^-o.oo 


0, oc:+2.43 

01.00_jg 7Q 




00+° °° 
'-'■'J'-' -0.00 


54, 


10+1-00 
■l'J-1.13 


54, 


Q5+0.58 

■^■-"-o-es 


Am. Py. 


0, 


.1 /im 


11 




2, 


70+1. 06 
■ '•^-0.91 


12, 


= 0+0-16 
■'"J-0.18 


00+° °° 




7 54+0.37 

'■^^^-0.37 




7, 


00+0.16 
,U»_o 17 


0, 


00+°°° 
■^"^-0.00 


Am. Py. 


2, 


,0 fim 


18 




21, 


70 + 1.01 
■ '■^-1.05 


15, 


04+0-42 
■'-'^-0.46 


3 50+°-92 




00+° °° 
'-'■'^'-'-0.00 




0, 


15+0.20 
■I'^-O.ll 


7, 


95+0.27 
■■^'-'-0.32 


Am. Py. 


5, 


,0 fini 


0, 


77+6.04 
•' '-0.75 


3, 


15+1.43 
.iO_j 27 


0, 


og + 0.00 

■'^^-o.oo 


qi+1.05 
■'•^^-1.57 




00+° °° 
'J-^'J-0.00 




0, 


OQ+O.OO 

■^'^-o.oo 


0, 


00+°°° 

■^*J-o.oo 


Fors. 


0, 


,1 /im 


1 


69+0 03 
•"^-0.04 


2, 


00+0. 04 
■■^°-0.04 


1, 


Q7+0.05 
-0.04 


1 10+°-'"' 
l-l'J-0.07 




1 14+0-04 
-^■-^^-0.04 




0, 


09+0.02 

■■^^-0.02 


0, 


gg + 0.03 

■"^-0.03 


Fors. 


2, 


,0 /im 


1, 


09+0.27 
■'-'^-0.39 


0, 


09+0-34 
■"^-0.09 


0, 


04+0-17 
■'-'*-0.04 


Q 07+0. 27 

U-U(_g g7 




00+° °° 
^■'^^-o.oo 




0, 


OQ+O.OO 

■^'^-o.oo 


0, 


34+°" 

"-•^-O.ll 


Enst. 


0, 


,1 /im 


1 


94+0.07 


2, 


qQ+0.08 
■■^"-0.07 


2, 


40+0. 05 
■^■^-0.06 


n 07+0. 18 

U-87_o gg 




1 54+0-09 
l-^*-0.06 




0, 


qo+0-03 
■•-'^-0.03 


0, 


34+0.04 
■'-'^-0.04 


Enst 2, 


,0 /im 


5, 


09+0.22 
■^^-0.23 


4, 


04+0-21 

■"^-0.24 


1, 


.7+0.22 
-0.18 


97+0.81 
' -0.38 




9 10+0.33 
^■l«-0.37 




0, 


OQ+O.OO 

■^'^-o.oo 


0, 


00+°°° 
■^*^-o.oo 


Silica 0, 


,1 /im 





Q1+0.03 
•"^-0.02 


0, 


89+0-03 
■°"-0.04 


0, 


49+0.02 
■^^-0.02 


O-SStaoi 




0.35t°°I 




0, 


QQ+0.00 

■^'J-o.oo 


0, 


70+0.02 
■ "J-0.02 


Silica 2, 


.0 /im 


2, 


9Q+0.13 
■^^-0.12 


4, 


09+0.22 
■°^-0.20 


3, 


c:7+0.13 
■^'-0.13 


0.85lJ:?° 




2.49t°:lJ 




0, 


OQ+O.OO 

■'-"^-0.00 


1, 


90+0.07 
■^"-0.07 


Silica 5. 


,0 iim 


9, 


40+0-34 


6, 


09+0.42 
■°^-0.61 


0, 


otr+0.46 
■OO_0.44 


0.00t°:~ 




10+0.44 
'J-i^-o.ia 




0, 


97+0.15 
■^'-0.19 


6, 


79+0.24 
■ '^-0.29 


PAH 6.2 


022+° °°2 


oq>;+0 001 


OQS+O-002 
U.098_o gg2 


369+° °°^ 





349+0-007 
■■'^''-0.007 


0, 


-,90+0.002 

•-^^^-0.002 


540+0-002 

U.O'lU_g 0Q2 


PAH 7.7 


000+" °°° 


004+°-°°^ 


064+° °°2 

U.U01_Q ggj 


900+0. 065 





OQo+0.000 
■'J^'J-0.000 





001 +0.002 

■'-'•-'^-0.003 


506+° °°^ 
U.aUD_g 023 


PAH 8.6 


069+° °°^ 


099+0 003 


101+°°°2 

U.iUl_Q gg2 


"S01+°159 





Q, -+0.034 
■'^^^-0.013 


0, 


Qog+0.000 
■'-"^'J- 0.000 


499+0-005 
u.'i»»_g 005 


PAH 11.3 


041+°-°°'' 


005+°-°°* 


059+° °°'' 

U.UOa_g ggj 


702+° °2'3 

U./UZ_g g3g 





QQ,+0.034 

■^"^i-o.ooi 





9 = 0+0.007 
•^"^"-0.006 


QQO+O.OOS 

u.ooo_g 006 


PAH 12.7 


000+° °°^ 

U.UUU_(j ggg 


000+°-°°° 

U.UUU_g QQQ 


044+0 003 

U.U*1_g ggg 


ot:o+0.014 





004+0.030 

■■-"^^-0.024 


0, 


994+0.003 
•^^^-0.003 


Q -,71+0.004 
'-'•-^ ' -^-0.004 



- 33 - 



Table 6. Same as TableU 



HD50138 HD58647 HD72106 HD85567 HD95881 HD98922 HD100546 









88.99 




15.01 


95.54 


42.40 




42.43 




46.95 




168.92 


Am. Ol. 


0. 


1 fim 




0, 




44 27+°- 
^^•^'-0.67 


10 11 +''■32 
J^U.ll_j^0.11 


0, 


00+°-°° 
■"'^-o.oo 


0, 


00+°-°° 
■'-"^-0.00 


0, 


09+0-78 
■"''-0.09 


Am. Ol. 


2, 


fim 


64.00l^'2g 


31, 




00+°-°° 


15 29+"-°^ 


70, 


Qg+0.84 
■^°-1.64 


0, 


00+°°° 
■'^"-o.oo 


82, 


0-1+0.51 
■■'^-0.76 


Am. Ol. 


5, 


,0 fim 


7 ig+13:75 
'•^^-6.18 


56, 


oc: + 1.53 
■■''^-1.79 


00+° °° 
'-'■uu_o 00 


c;9 70 + 39.11 


0, 


gq+5.57 
■"^-0.69 


89, 


70+0.44 

•'■^-0.50 


0, 


-,Q+1.28 
■^^-0.19 


Am. Py. 


0, 


.1 /im 


00+" °° 
U-UU_o 00 


0, 


00+°°° 

■'^^-o.oo 


10 79+°-44 


77+1.06 
U- ' ' -0.77 


0, 


nq+0.37 
■"■^-0.03 


0, 


00+°-°° 
•^'^-o.oo 


0, 


00+°°° 

■^*^-o.oo 


Am. Py. 


2, 


,0 fim 




0, 


00+°°° 
■'^'^-0.00 


6 45+°-5i 


5 qq+2.09 


0, 


ic:+0.24 
■^°-0.12 


0, 


00+°-°° 
•^'J-o.oo 


0, 


00+°°° 
■'^'^-o.oo 


Am. Py. 


5, 


,0 fini 


27+1-°^ 


0, 


00+°°° 
■'^'^-o.oo 


00+°-°° 


00+° °° 

n 03+0.35 

u.»o„o.49 


0, 


00+°°° 
■^^-0.00 


0, 


00+°-°° 
•^'^-o.oo 


0, 


00+°°° 
■^*J-o.oo 


Fors. 


0, 


,1 /im 


-1 99+0.11 
^•^^-0.11 


0, 


64+004 
■"^-0.04 


7 67+°-i° 
'■°'-o.io 


1, 


47+0.09 
-0.08 


0, 


c:o+0.04 
•^''-0.04 


5, 


go+0.14 
■"'-'-0.13 


Fors. 


2, 


,0 /im 


oo+°-i® 

^^•'JU_o 00 


0, 


00+°°° 

■'^'^-o.oo 


00+° °° 
>J-UU_o 00 


0.001°:°^ 


0, 


00+°°° 
■^'-'-0.00 


0, 


00+°-°° 
•^'^-o.oo 


0, 


00+°" 

■'^^^-o.oo 


Enst. 


0, 


,1 /im 


1 49+0.17 
^•^^-0.20 


0, 


99+0.06 
■^^-0.06 


6 74+° " 
'^-0.17 


1.68ir,i 


2, 


70+0.15 

■ '°-0.12 


0, 


Qq+0.07 
•^■'-0.07 


0, 


91+0.12 
■^^-0.12 


Enst 2, 


,0 /im 


o qo + 0.59 
»-3»_0.88 


1, 


Qq+0.19 
"■^-0.15 


i7.i8t°:^^ 


3 46+^ °'' 


7, 


C7+0.28 
-0.27 


5, 


gg+0.18 
•°^-0.15 


4, 


00+0. 29 
■■^■^-0.27 


Silica 0, 


,1 /im 


^•00-0.32 


0, 


7-1 +0.03 
■ '^-0.03 


5.92i°-°«9 


0.82+°-^6 


3, 


90+0.08 
■^=^-0.08 


0, 


07+0.03 
•■^'-0.02 


0, 


00+0.07 
■*'°-0.07 


Silica 2, 


.0 /im 


^ no + 0.43 


0, 


00+°°° 
■'^'^-0.00 


0.98t°:}« 


2 56+1-1* 


3, 


y,c:+0.35 
■^■^-0.31 


0, 


gQ + 0.02 

•^'^-o.oo 


6, 


00+0. 14 
■■'°-0.14 


Silica 5. 


,0 iim 


Q QQ+l-SS 

».03_2 20 


8, 


4,5+0.33 

■^■^-0.27 


00+° °° 


5 61+* °'^ 
0-"l-5.61 


10, 


40+0.94 
■^■^-0.87 


2, 


gg + 0.37 

•"^-0.33 


0, 


00+°°° 

■'^^-o.oo 


PAH 6.2 




061+° °°* 
U.UDl_o 003 


„ t:-, 1+0.002 
u.oii_o 002 


n c9c:+0.026 
u.ozo_o 020 


1 264+°-°l° 
^■^"^-0.010 










PAH 7.7 


603+°-3^^ 


n -1 7Y+0.003 
'^•^ ' ' -0.003 


n c;-, + 0.002 
u.olo_o 003 


n K77+0.021 
'-0.029 


1 K9S+0.006 
1.0ZO_o 005 


1 505+0.091 
i.oao_o 091 


71 1 +0.255 

' ^^-0.266 


PAH 8.6 


5 oo';+°"i° 

O.UUO_o 


U.lSO_o 004 


r, 07-1+0.002 
'J-'Jl J^_0.002 


n /1S9+0101 

U-48i_o 139 


l-346l°:»^ 


9 9of;+0.049 
2-^»O_0.045 


5fiQ+0140 
0. JU»_Q 124 


PAH 11.3 


Q o'3/i+0.213 
0.00'i_o 205 


U.089_o 006 


r, 900+0.007 
U.zoo_o.o08 


191+0 024 
u.l»l_Q 026 


0.863t°:°r7 


2 030+° °*° 
z.uou_Q 056 


4 116+°-^*'' 
'l.liD_0 3g4 


PAH 12.7 




063+° °°=' 
U.UDO_o 003 


n ns'^+0-002 

O.U»3_o 002 


150+°"2 
U.10U_o 096 


0.286l°:»^ 


1 654+°-°5'' 
i.DO'i_o 070 


2 369+° "'' 
z.oo»_Q 140 
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Table 7. Same as TableU 



HD101412 HD104237 SS73 44 HD139614 HD142666 HD142527 HD144432 











62.55 




117.56 


41.20 


34.40 


19.94 




210.91 




77.59 


Am. Ol. 


0. 


1 fira 


0. 


00+"°° 

'^'^-0.00 


0, 


00+°°° 


69+°-l2 


Q 17+0,32 

"■^'-0.14 


7 99+0.97 
' •^^-0,93 


37, 


14+0,70 

-^-^-0.62 


5, 


Q9+3.26 
■='^-2.89 


Am. Ol. 


2, 


/im 


63, 


^^5-3.35 


0, 


00+°-°° 
■^^-0.00 


00+° °° 
"■""-0.00 


38.52_-^ Qg 


45.20t^:^* 


13, 


c:7+0,57 
""-0.52 


42, 


oo~l~2 . 32 
■'-"'-2,61 


Am. Ol. 


5, 


,0 lira 


4, 


on+4.32 


73, 


co+0.82 
■"-'■^-0.79 


Ql 45+0-23 


C7 00+0,86 


49+3.08 
•-"•^^-3.08 


0, 


00+°°° 
■""-0.00 


0, 


Q1+0.07 
■"^-0,01 


Am. Py. 


0, 


.1 lira 


0, 


QQ+0.00 

■*^^-o.oo 


0, 


00+°°° 
■'^'^-0.00 


r, 01 +0.06 

^■■'^-o.oe 


00+° °° 
"•""-0.00 


00+° °° 
"■""-0,00 


0, 


00+°-°° 
■""-0.00 


19, 


04+1-01 
■"^-1.05 


Am. Py. 


2, 


,0 ixra 


0, 


00+°°° 
■'-"^-0.00 


10, 


70+1-86 
■ "-'-1.79 


4 77+0.18 
*■ ' ' -0.20 


00+° °° 
"•""-0.00 


c QC + 1,41 

6.36_j 47 


22, 


gg + 0.51 

■^°-0.49 


14, 


Qtr+0.59 
■'^'•'-0,63 


Am. Py. 


5, 


,0 /im 


0, 


00+°°° 
■^'^-o.oo 


3, 


Q9+2.48 
■^^-2.68 


00+° °° 
"■""-0.00 


00+° °° 
"•""-0.00 


22 15+2-™ 
zz.iO_2 40 


0, 


00+°-°° 
■""-0.00 


0, 


QQ + O.IO 

■""-0,00 


Fors. 


0, 


,1 ura 


2, 


45+0.13 

■^'^-0.13 


2, 


QQ+0.16 

■'J'^-0.17 


Q 00+0.01 

"■'-'•-'-0.02 


90+° °" 
"•^"-0.04 


1 02+°-°** 
-'-■"^-o.os 


5, 


7n+0.07 
■'"-0.09 


2, 


qfi+0.04 
■■^"-0,05 


Fors. 


2, 


,0 /^m 


0, 


00+°-°° 
■'-"^-0.00 


0, 


00+°°° 
■""-0.00 


00+° °° 
"■""-0.00 


n -is+0,22 

"•1^-0,13 


00+° °° 
"■""-0,00 


0, 


00+°-°° 
■""-0.00 


0, 


00+°°° 
■""-0,00 


Enst. 


0, 


,1 /^m 


1, 


0-1+0.15 

■■^^-0.18 


1, 


1 9+0.22 
■^'^-0.23 


45+0 02 
"■^°-0.02 


50+° °" 


1 00+°-" 
^■""-0,13 


3, 


90+0.15 

■^''-0.12 


1, 


09+0.08 

■'^^-0,08 


Enst 2, 


,0 /im 


20, 


on+0-56 


6, 


40+0-57 
■^"-0.41 


00+°-°° 
"■""-0.00 


-I 79+0,19 
^•'^-0,20 




6, 


97+0.20 
■^'-0.22 


4, 


0-1+0.15 


Silica 0, 


,1 /^m 


3, 


7fi+0.13 
■'°-0.18 


1, 


00+0. 09 
■■^"-0.06 


00+°-°° 
"■""-0.00 


n 00+0,02 

"•■-''-'-0,01 


o.o8«:°^ 


6, 


trfi+0.06 
■''°-0.05 


0, 


in+o.05 

■^"-0,05 


Silica 2, 


.0 lira 


3, 


00+0.20 
■°'^-0.22 


0, 


50+0.15 
■^"-0.17 


00+°-°° 
"■""-0.00 


29+° °^ 
"•^^-0,05 


00+°°° 
"■""-0,00 


4, 


c:c; + 0.13 

■°^-o.io 


1, 


Q7+0.28 
■=" -0,28 


Silica 5. 


,0 /im 


0, 


00+°°° 

■^'^-o.oo 


0, 


47+1-28 
-0.47 


00+°^°° 
"■""-0.00 


06+°-2^ 
"•"°-0,06 


12 46+^ °'* 


0, 


00+°°° 
■""-0.00 


7, 


1tr+0.61 
■^^-0.61 


PAH 6.2 


Ol C-+0.005 


Q -,02+0.009 

^•^^^-0.011 


007+° °°^ 
"■""'-0.001 


504+° °°3 

U.OU^_Q 002 


565+0 006 


1 Qn9+0-005 
l.ouz_o 006 


206+° °°^ 

U.ZUD_o 004 


PAH 7.7 


o 077 + 0.002 

' -0.002 


o.ooot°:°°° 


001+° °°^ 
"■""^-0.001 


464+° °°2 
"•*°^-0.002 


„ 5-,, +0.007 
U.01i_o 008 


1 as-^+0-006 
l.ooo_o 006 


0.000t°:°°° 


PAH 8.6 


oqc: + 0.005 

U.OOO_g QQg 


0.863t°:°f, 


000+° °°^ 


oQl +0-005 


000 + 0.011 

u.ooo_o 012 


1 070+0-012 
l.»lu_o 010 


0.048t°-°3f 


PAH 11.3 


n irrn + 0.005 


0.602+°-}™ 


000+° °°° 
"■"""-0.000 


(, 051+0.013 

U.001_o 010 


502+° °^° 
U.OU/_o 017 


Q47+0.035 
"•"^' -0,032 


169+°°^^ 

U.lDM_o 024 


PAH 12.7 






n ns'^+O-002 
U-U»3_o 002 


260+° °°'^ 


Q 1 -,7+0.013 
"■^^'-0.010 


n 109+0. 008 

u.loz_o 007 


042+° °l2 
u.u^z_o 012 
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Table 8. Same as TableU 



HD144668 HD145263 HD150193 HD152404 51 Oph HD163296 VV Scr 









65.92 




58.18 


221.96 




38.74 




65.46 




80.59 


35.88 


Am. Ol. 


0. 


1 fira 


00+" °° 




7 51+1.57 


0.011°-? 


53, 


7fi+3,03 
'"-2.69 


0, 


00+°-°° 
■"'^-o.oo 


0, 


00+0,00 
■""-0.00 


00+° °° 
"•""-0.00 


Am. Ol. 


2, 


/im 


71 Qo+2.73 
1 ±.oo_3 4y 


58.64t^; 


56.5lll:f8 


11, 


77+2.26 
' ' -2.45 


45, 


qq+0.00 

■■'■5-0,00 


62, 


3g+l,81 
■■'"-2.30 


21.90+i°3f 


Am. Ol. 


5, 


lira 


9 1 9+8.57 
1 S8 




- oc:+4.26 
'■•53-4.43 


'^■^^-0.09 


0, 


Qi+0,12 
■^1^-0.01 


0, 


00+°°° 
■'^'-'-0.00 


1, 


qi+5.33 
■■5^-1.25 


38.89tl°:« 


Am. Py. 


0, 


.1 lira 


00+°'°° 
u.uu_o 00 




00+° °° 


00+° °° 
u.uu_o 00 


21, 


74+1-11 
■'^-1.20 


0, 


00+°°° 
■^^-0,00 


1, 


C7+O.39 
■■"-0.36 


0.001°:°° 


Am. Py. 


2, 


,0 lira 


02+1°^ 




Q 04+102 


20 76+°-™ 


0, 


t;9+0.72 
■"^^-0.43 


0, 


00+°°° 
■'^'-'-0,00 


18, 


44+0.77 
■*^-0.90 


0.00lg:°° 


Am. Py. 


5, 


,0 lira 


02+°'^^ 




OQ+°-l^ 


00+°-°° 


0, 


00+°°° 
■'-"^-0.00 


33, 


60+°°° 
■"'-'-0,00 


0, 


00+0.00 
■""-0.00 


24.5lll:« 


Fors. 


0, 


,1 lira 


2 qi+°13 
^•»J^_0.17 




r, 71+0.13 

^- ' ^-0.13 


9 ss+0.05 
^■*'**-0.05 


2, 


11+0.05 
■^^-0.06 


0, 


9q+0,00 
■^■5-0,00 


2, 


gi+0.07 
■"^-0.10 


1 10+0.21 

^■1^-0, 15 


Fors. 


2, 


,0 lira 


'^•'^^-0.02 




01+°-^" 


2 46+°-^° 
^■™-0.33 


0, 


00+°°° 
■'-"^-0.00 


3, 


19+0,00 

■^^-0,00 


0, 


ni+o,20 

■"1-0.01 


27+0-37 
"■^'-0,25 


Enst. 


0, 


,1 lira 


"•"^-0.22 




1 05+°-25 
l.UO_o 26 


1 79+°- 1° 


1, 


qn+o-os 
■■^'^-o.os 


0, 


90+0,00 

■^"-0,00 


1, 


OP-+0.14 
■"3-0.13 


63+° °'' 
"■"■J-0,08 


Enst 2, 


,0 lira 






5 4Q+0.54 


q 7Q+0.28 
■5- '='-0.30 


4, 


qg+0.21 
■"''-0.22 


1, 


7Q+0,00 
■ '"-0,00 


2, 


t:c; + 0.36 
■"3-0.34 


q 7q + 0.90 
■'■ '■'-0.65 


Silica 0, 


,1 /^m 


2.08l°:l^ 




64+°-^* 


1 07+0.05 
-"-■^'^-0.05 


0, 


-,1+0.07 
■-^-^-0.06 


0, 


■"^-0,00 


0, 


f,o+0.06 

■"^'-o.oe 


u.Da_o 21 


Silica 2, 


.0 lira 


2.59l°jl 




- c:7+0.40 
'■^'-0.46 


4 96+°-2l 


3, 


gg+0.18 
■°°-0.20 


0, 


00+°°° 
■"'^-0,00 


4, 


c:c; + 0,27 
■33-0.33 


00+° °^ 
"■""-0.00 


Silica 5. 


,0 /im 


10.97irj 




00+° °° 


4.»»_o.g4 


0, 


09+0-12 
■^^-0.02 


14, 


7-, +0,00 
■ '^-0,00 


4, 


00+1.08 
■""-1.13 


''■l"-2.42 


PAH 6.2 




0, 


006+°°°2 
■'-"^°-0.002 




n r|7c:+0.001 
U.U(0_o 002 


000+° °°° 
u.uuu_o 000 






041+0.004 
U.O'il_o 004 


PAH 7.7 


065+° °'^'' 





Q26+0.001 
•'-'^"-0.001 


000+° °°° 
u.uuu_o 000 


000+° °°° 

U.UUU_o 000 


012+°-°°° 
u.uiz_o 000 


U.lzu_o 100 


307+0.016 
"•■'"'-0.012 


PAH 8.6 


640+°-°'"' 

U.D'iU_o 048 


0, 


ntr9+0.002 
•'^'^^-0.002 


1 664+°-°'^l 
l.DD1_o 047 


030+° °°* 
U.UOU_o 008 


n /iqo+0,000 
u.'ioo_o 000 


1 1 30+0.092 
l.lou_o 059 


n qc:K+0.044 
"••330_o 039 


PAH 11.3 


792+°-°43 
u. oz_o 048 





077+0.002 
-0.002 


027+° °='^ 
'-'■U^'_0.023 


010+° °°^ 
U.UiU_o 005 


n fi7s+o.ooo 

U.bC»_o 000 


591+° °™ 
u.j»±_o 070 


147+0.026 
"•1^^' -0.020 


PAH 12.7 


r, 171 +0.062 
'^■^'^-0.057 


0, 


QOO+O.OOl 
•'J '^'-'-0.000 


r, qr,r, + 0.031 

u.ozz_o 037 


043+° °°3 
U.U'iO_o 003 


169+°°°° 
u.iua_o 000 


ig7+0.039 
U-ID' -0,036 


0.058_o;oi7 



- 36 - 



Table 9. Same as TableU 



T CrA 



HD179218 



WW Vul 



HD190073 



HD203024 











23.73 




109.27 




50.51 




100.87 


254.32 


Am. Ol. 


0. 


1 fim 


31. 




31, 


0-+1.41 
,»<'_1 72 


0, 


00+0.00 
■'^'-'-0.00 


0, 


00+0.00 


5 50+49,75 

'^■'''^-5.53 


Am. Ol. 


2, 


fim 


54, 




32, 


■"^-1,57 


50, 


50+0.29 
""^-0.28 


31, 


,g+2,27 
1°-1,71 


95 01+3.06 
ZO.Ui_j^y 32 


Am. Ol. 


5, 


,0 /im 


0, 


r, -I +0.54 

^^-0.17 


0, 


on+0.64 
■•^^-0.24 


0, 


00+0.00 
■'^'-'-0.00 


20, 


90+2.99 
■^^-3.66 


46 25+^-^'' 


Am. Py. 


0, 


.1 /im 


0, 


00+°-°° 
■'^'^-0.00 


0, 


00+°°° 
■'^'^-0.00 


18, 


09+0.15 
■'^''-0.16 


0, 


00+0.00 

■^'^-o.oo 


7 qQ+12.54 

' -^^-1.39 


Am. Py. 


2, 


,0 /im 


10, 




2, 


Qc+0,60 
■^■^-0.63 


13, 


99 + 0.18 
■^^-0.21 


0, 


00+0.00 
■'-"^-0,00 


6.85t°i^ 


Am. Py. 


5, 


.0 /im 


0, 


00+°°° 

■^'^-o.oo 


0, 


00+°°° 

■'^^-o.oo 


0, 


ni +0.08 
■'^l-O.Ol 


16, 


-,g+l,02 

■-'-"-0,87 
Og+0,05 
■'-"^-0,05 


0.00t°:°° 


Fors. 


0, 


,1 /im 


0, 


■o'-o.io 


0, 


90+0 09 
■'''-'-0.09 


1, 


4,5+0.03 
■^'-'-0.03 


1, 


1 98+^-*^ 
1-^^-0.16 


Fors. 


2, 


,0 /im 


0, 


00+°-°° 
■'-"^-0.00 


0, 


00+0.00 
■'^'^-0.00 


0, 


46+0.13 
■^"^-0.14 


1, 


90+0,22 
■^■^-0,18 


00+° °° 
'-'■'^^-0.00 


Enst. 


0, 


,1 /im 


0, 


t.7+0,20 

■^'-0,22 


6, 


00+0.19 
■'^"-0.18 


1, 


24+0.04 
■^^-0.04 


2, 


06+°°° 
■'-"^-0,08 


67+^-^^ 
'-'■'^'-0.18 


Enst 2, 


,0 /im 


0, 


■"'^-0.36 


17, 


oc+0.26 
■■''^-0.30 


3, 


57+0.11 
■0' -0.16 


10, 


09+0.35 
■'^^-0,28 


4 83+6.40 
^■^■^-0.71 


Silica 0, 


,1 /im 


0, 


trn + 0,08 
■0^-0,09 


6, 


00+0,10 
■^O-O.ll 


0, 


09+0.02 
■'^^-0.02 


1, 


7^+0,08 
■'^-0,07 


1 97+4,28 
-'-■^' -0.48 


Silica 2, 


.0 /im 


2, 


-,0 + 0.21 
■^"-0.24 


0, 


00+0.02 
■'^'^-0.00 


0, 


oc+O.lO 
■''°-0.09 


0, 


00+0.00 

■'J'^-0,00 


n 90+2.06 

'-'■^■-'-0.23 


Silica 5. 


,0 iim 


0, 


01+0.54 
■'^^-0,01 


0, 


00+0.00 
■'J'J-0.00 


10, 


gO+0.30 


16, 


9n+o.88 

■^'^-0.72 


00+° °° 


PAH 6.2 










n5i+° °°i 


,-, ,97+0.003 

'^•^^'-0.003 


029+° °°^ 

U.UZ»_g QQ7 


PAH 7.7 


000+" °°° 


3.286t°;}f, 


000+° °°° 

U.UUU_p QQQ 


o.ooot°°S° 


002+°-°^5 

U.UUZ_g QQ2 


PAH 8.6 


090+° °^^ 


q oc:q+0.090 

0.003_Q ggg 


000+° °°° 

U.UUU_Q QQQ 


0.286t°°}i 


0.028l°:2g 


PAH 11.3 


250+°-°32 


9 sc;7+0-074 


055+° °°^ 

U.UOO_Q QQ3 


432+° °^* 




PAH 12.7 


016+°-°l'^ 

U.UlD_(j 


+0.036 
^■1 '3-0.036 


000+° °°^ 

U.UUU_g QQQ 


092+° °°* 
u.uaz_Q 


o.288inr7 
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Table 10. Best fit dust parameters and the reduced of tfie fit (17-35 /um). Listed are the 
derived mass fractions of each dust component in %. 









AB Aur 


HD31648 


HD35187 


HD35929 


HD36112 


HD244604 


HD36917 










58.93 




46.29 




31.46 




27.89 




49.36 


37.61 




187.87 


Am. Ol. 


0. 


1 fim 


30. 




0, 


00+°-°° 
■'^'^-o.oo 


0, 


24+0.69 
■^^-0.22 


0. 


00+°°° 
■'-"^-0.00 


0, 


45+3.62 
■^^-0.45 


00+°-°° 
'^•'-"J-0.00 


0, 


00+°°° 
■'J^-o.oo 


Am. Ol. 


2, 


,0 fim 


0, 


00+°°° 


0, 


00+°°° 
■'^^-o.oo 


0, 


n^+1.20 


14. 


-,0+2.72 
■^^^-2.84 


0, 


00+°°° 
■^'-'-0.00 


00+° °° 
^•*J'J-o.oo 


0, 


00+°°° 
■'^^-o.oo 


Am. Ol. 


5, 


,0 /im 


0, 


00+°-°° 


0, 


QQ+4.90 
■^■^-0.93 


58, 


43+1-34 
■^■-'-1.64 


9. 


11+1.17 
■^^-3.52 


0, 


00+°°° 
■'^'-'-0.00 


00+° °° 


0, 


00+°°° 
■'^'-'-0.00 


Am. Py. 


0, 


.1 /im 


64, 


■«'-1.28 


84. 


0-+3.56 
"^-3.03 


1, 


07+3. 50 
■°'-1.43 


71, 


gg + 2.98 

■"''-3.23 


90, 


4g+0.44 
■*°-1.03 


84.6212-53 


97, 


ie+0.30 
■^"-0.30 


Am. Py. 


2, 


,0 /im 


0, 


00+°-°° 


1, 


cq+6.51 
■"■^-1.63 


37, 


17+1.51 
' -2.09 


1, 


,1+3.56 
■^^-1.06 


0, 


00+°°° 
■'^'-'-0.00 


"■^■-'-0.23 


0, 


00+°°° 

■'J^-o.oo 


Am. Py. 


5, 


.0 iim 


0, 


00+°°° 


6, 


Qg+3.46 
■'-'"-3.90 


0, 


14+0.33 
■-^^-0.11 


0, 


00+°°° 

■'^'^-o.oo 


0, 


QQ+0.02 

■'^^-o.oo 


c: Qc: + 3.36 
3.oO_2.43 


0, 


00+°°° 
■'^'^-o.oo 


Fors. 


0, 


,1 /im 


0, 


07+0.08 
■^'-0.05 


2, 


OC-+0.12 


0, 


00+°°° 

■'J^-o.oo 


0, 


qo+0.12 
■■^■^-0.14 


1, 


p;n+0.07 
■^^-0.07 


4 17+0.12 
^•I'-O.ll 


0, 


00+°°° 
■'^'-'-0.00 


Fors. 


2, 


.0 /im 


0, 


r,o + 0.08 

■^"-0.10 


0, 


gq+0.10 
■"^-0.10 


0, 


00+°°° 
■'^'J-o.oo 


0, 


1 7+0.23 
■^'-0.14 


2, 


QQ+0.12 

■^^^-0.13 


00+° °° 
'^•'-'^-0.00 


2, 


04+0.29 

■^4_o.3l 


Enst. 


0, 


.1 jim 


0, 


00+°-°° 
■'-"^-0.00 


1, 


c:n+0.09 
■^'-'-0.08 


0, 


no+o.io 

■^^^-0.02 


3, 


1Q+0.40 
■^"-0.36 


1, 


og+0.11 
■■^"-0.09 


9 44+0-15 
^•^^-0.17 


0, 


00+°°° 

■'J^-o.oo 


Enst 2, 


.0 iim 


1, 


^1+0.11 
"^^-0.09 


0, 


00+°°° 
■^'-'-0.00 


1, 


41+0.16 

■^-^-0.13 


0, 


Q1+0.27 

■"J^-o.oi 


0, 


00+°°* 

■'^^-o.oo 


00+° °° 


0, 


00+°°° 
■'^^-o.oo 


Silica 0, 


,1 /im 


0. 


00+°°° 
■'^'^-0.00 


1, 


■■'"-o.oe 


0, 


oQ+0.07 
■"•-'-0.08 


0, 


00+°°° 
■'^'J-o.oo 


0, 


OQ+0.05 
■^'-'-0.00 


2 62+°-" 
^•"^-0.19 


0, 


00+°°* 
■'^'^-0.00 


Silica 2, 


,0 /^m 


3, 


r, 7+0. 08 
■^'-0.10 


0, 


00+°°° 
■^'J-o.oo 


0, 


00+°-°° 
■^*J-o.oo 


0, 


00+°°° 
■'-'^-0.00 


2, 


07+0.11 

■^' -0.07 


06+°-2'' 
^•*J°-0.06 


0, 


00+°°° 
■'^^-o.oo 


Silica 5. 


,0 /im 


0, 


00+°°° 
■'-"J-0.00 


0, 


00+°°° 
■'^'-'-0.00 


0, 


94+0.35 
■^^-0.21 


0, 


00+°°° 
■'-"^-0.00 


1, 


90+0.35 
■^°-0.71 


00+° °° 

'J-^'J-o.oo 


0, 


00+°°° 
■'^'-'-0.00 



Table 11. Same as TablefTOl 



HD37258 BF Ori HD37357 HD37806 HD38120 HD250550 V590 Mon 











21.03 




38.37 




10.63 




74.21 




38.11 




16.20 




12.87 


Am. 01. 


0, 


1 fim 


0, 


0(1+0.00 
'-"J-0.00 


0, 


00+0.00 
■"'^-o.oo 


0, 


00+0.00 
■'J'J-o.oo 


5. 


7C+0.87 
■'°-1.02 


28, 


04+1.07 
■'-'^-1.11 


6. 


Qq+5.34 
■"■'-4.93 


0, 


00+0.00 
■'^"-0.00 


Am. 01. 


2, 


,0 /im 


0, 


00+0.00 
'-"^-0.00 


0, 


00+0.00 
■^'-'-0.00 


0, 


00+0.00 
■^"^-0.00 


0. 


00+0.00 

■'J'^-o.oo 


0, 


00+0.00 
■'^^-o.oo 


0. 


90+4.33 
■^"-0.28 


0, 


c.o+2.98 

■"S-o.ee 


Am. 01. 


5, 


,0 fim 


0, 


09+1.02 
■^^-0.02 


0, 


00+0.00 
■'^'-'-0.00 


1, 


59+2.18 
■^"-1.40 


0. 


00+0.00 
■'-"^-0.00 


0, 


97+0.42 
■^'-0.22 


19. 


7q+5.19 
■'■^-5.85 


21, 


47+4.14 

■^' -4.30 


Am. Py. 


0, 


,1 fim 


0, 


nK+2.63 
■*J'^-0.05 


93, 


trq+0.28 
■'^'-'-0.30 


64, 


g^+4.91 

■^^-5.32 


89. 


49+1.03 
■^^-0.87 


30, 


nQ+0.46 
■'-''-'-0.50 


28. 


gO+5.16 
■°"-5.82 


51, 


c:q+4.76 
■'-'■-'-4.22 


Am. Py. 


2, 


,0 fim 


95, 


01 +0.37 
■"^-0.87 


0, 


00+0.00 
■'^'-'-0.00 


26, 


05+5.48 
■■^^-5.71 


0. 


00+0.00 

■'^'J-o.oo 


0, 


00+0.00 

■'J'-'-o.oo 


30. 


9Q+5.36 
■^"-4.21 


18, 


Oq+6.13 

■"'-'-5.43 


Am. Py. 


5, 


,0 /im 


0, 


20+1-12 
■^"-0.18 


0, 


01+0.12 
■^1-0.01 


2, 


7t:+3.16 
■ '^-2.58 


0. 


-ic;+0.16 
■1^-0.10 


40, 


00+0.99 
■'-"'-0.99 


5. 


4q+9.30 
■^■'-5.01 


0, 


00+0.02 
■""-0.00 


Fors. 


0, 


,1 /im 


0, 


09+0.07 
■°^-0.08 


2, 


00+0.15 


0, 


tr7+0.09 
■^' -0.08 


3. 


90+0.08 
■^■^-0.08 


0, 


00+0.00 
■'"'-0.00 


0. 


00+0.00 
■"'^-0.00 


2, 


55+0.13 
■^°-0.12 


Fors. 


2, 


,0 /im 


0, 


91+0.15 
■^^-0.13 


0, 


■"^-0.17 


1, 


90+0.14 
■^"-0.13 


0. 


70+0.13 
■'^'-O.ll 


0, 


00+0.00 

■""-o.oo 


1. 


99+0.12 
■^^-0.10 


0, 


04+0. 18 
■^4-0.18 


Enst. 


0, 


,1 /im 


1, 


90+0.44 
■^■^-0.45 


0, 


10+0.20 
■^«-0.12 


0, 


47+0.30 
■^' -0.25 


0. 


47+0.10 
■*'-0.09 


0, 


84+0-07 
■"^-0.07 


0. 


ns+0-20 

■'J*'-0.07 


0, 


-1-1+0.28 
■11^-0.10 


Enst 2. 


,0 /im 


1, 


1 1+0.53 
■1^-0.51 


0, 


00+0.00 

■'^^-o.oo 


0, 


07+0.28 
■"' -0.39 


0. 


00+0.10 
■'^'^-o.oo 


0, 


00+0.00 
■"'-'-0.00 


0. 


cs+0.20 
■°"-0.22 


1, 


94+0.23 
■^^-0.37 


Silica 0, 


,1 /im 


0, 


4g+0.12 
■™-0.12 


3, 


cg+0.44 
■^°-0.45 


1, 


■'^"-0.14 


0. 


10+0. 04 
■^**-0.05 


0, 


oq+0.03 
■"■-'-0.03 


0. 


00+0.00 
■"'^-0.00 


2, 


79+0.16 
■ '^-0.16 


Silica 2, 


,0 /ini 


0, 


00+0.00 
■'-'^-0.00 


0, 


00+0.00 
■^^-0.00 


0, 


00+0.00 

■'^^-o.oo 


0. 


00+0.11 
■'J'^-o.oo 


0, 


00+0.00 
■""-0.00 


2. 




0, 


nq+0.43 
■"'-'-0.03 


Silica 5, 


,0 iim 


0, 


Oc+0.35 
■'-'^-0.04 


0, 


00+0.00 

■'^'-'-0.00 


0, 


1 9+0.60 
■1^-0.11 


0. 


00+0.00 
■'-"^-0.00 


0, 


00+0.00 
■'^"-0.00 


4. 




0, 


00+0.00 
■""-0.00 
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Table 12. Same as TablefTOl 



HD50138 


HD58647 


HD72106 


HD85567 


HD95881 


HD98922 


HD100546 











106.34 


38.52 




42.33 




78.85 




152.74 




180.41 




67.64 


Am. Ol. 


0. 


1 fim 


37, 


1^+0.89 


o (,r5+22.19 


8, 


07+1.84 


44, 


27+0,73 
'°'-0.67 


0, 


(,(,+0.03 
■"'-'-0,00 


52, 


7-1 +3,39 
' J^-1.74 


0, 


(,(,+0.00 
■'^"-0.00 


Am. Ol. 


2, 


fim 


0, 


'-"^-0.00 


91 sq+''-^* 

^i-Oa_g 00 


0, 


(,(,+0.00 

•*J'J-o.oo 


0, 


(,(,+0.00 
■'J'J-0.00 


0, 


07+0.73 
■'-" -0.07 


0, 


(,(,+0.00 
■'J'J-0.00 


0, 


(,(,+0.00 
■""-0.00 


Am. Ol. 


5, 


,0 /im 


1, 


(5(5+0.77 
"0-0.81 




7, 


00+1.02 
•■^°-1.25 


0, 


(,(,+0.00 
■'^'J-0.00 


0, 


(,(,+0.00 
■'^'-'-0.00 


0, 


■oolHS 


0, 


(,(,+0.00 
■""-0.00 


Am. Py. 


0, 


.1 /im 


52, 


(57+0.47 
■O'-0.51 




79, 


Qq+0.81 
■^■^-0.84 


48, 


(,p.+0.65 
■'^'^-0.70 


88, 


7c:+0.15 
■ '"-0.28 


41, 


■o&l 


91, 


70+0.17 

■'°-0.24 


Am. Py. 


2, 


,0 /im 


0, 


(,(,+0.00 


Ol -, 7+8.82 
^-'-■-'-'-7.83 


0, 


(,(,+0.00 
•'-"J-0.00 


0, 


(,(,+0.00 
■'-"J-0.00 


0, 


(,(,+0.00 
■'^'-'-0.00 


0, 


(,(,+0.00 

■"'J-0.00 


0, 


(,(,+0.00 
■""-0.00 


Am. Py. 


5, 


.0 /im 


0, 


09+0.18 
■"^-0.02 


'-'■'^*J-o.oo 


0, 


(,(,+0.00 

•*J^-o.oo 


0, 


(,(,+0.00 
■'-"^-0.00 


0, 


(,(,+0.00 
■'^'J-0.00 


0, 


(,(,+0.00 
■"'J-0.00 


0, 


(,4+0.25 
■'J*-0.03 


Fors. 


0, 


,1 /im 


2, 


-,(,+0.04 


n s3+°-i2 


1, 


Qc+0.04 


1, 


-, 7+0.06 
■^'-0.05 


3, 


40+0.05 

■^■^-0.07 


1, 


-,7+0.04 
■-^-'-0.07 


4, 


r,-,+0.15 

■^J^-0.11 


Fors. 


2, 


,0 /im 


1, 


qi +0.08 
■■^J^-0.09 


o 4Q + 0.21 
■^■^^-0.22 


0, 


(,(,+0.00 

•*J'J-o.oo 


0, 


5(,+0.11 


0, 


(,^+0.10 
■'-'"-0.04 


0, 


(,(,+0.00 
■'J'J-0.00 


0, 


Q-,+0.14 
■^J^-0.16 


Enst. 


0, 


,1 /im 


2, 


-,(5+0.15 


^■'J^-0.49 


1, 


Q7+0.11 
•*''-0.11 


3, 


qtr + 0.06 
■■^3-0.09 


7, 


(,(5+0.09 
■'^"-0.11 


4, 


(54+0.05 
■"^-0.06 


0, 


(,(,+0.00 
■""-0.00 


Enst 2, 


,0 /im 


1, 


9n+0-20 
■^^-0.20 


n Q1+0-61 
'-'■^-'--0.54 


0, 


(,(,+0.00 

•*J^-o.oo 


0, 


(,(,+0.00 
■""-0.00 
(,(5+0.05 
■'-'"-0.03 


0, 


(,(,+0.00 
■""-0.00 


0, 


(,(,+0.00 
■'J '^-0.00 


0, 


(,(,+0.00 

■'^"-o.oo 


Silica 0, 


,1 /^m 


1, 


74+0.05 


9 0-1+0.24 
^■•^-'--0.29 


0, 


(,(,+0.00 
•'-"J-0.00 


2, 


0, 


(54+0.08 
■"^-0.07 


0, 


4-1+0.03 
■*4-0.03 


3, 


(,0+0.11 

■'^"-0.10 


Silica 2, 


.0 /jm 


0, 


(,(,+0.00 
■'-"^-0.00 


"■""-0.00 


0, 


(,(,+0.00 
•'J'J-0.00 


0, 


(,(,+0.00 

■'JU-o.oo 


0, 


(,(,+0.00 

■""-0.00 


0, 


(,(,+0.00 
■'J'J-0.00 


0, 


09+0.28 
■'-'^-0.02 


Silica 5. 


,0 iim 


0, 


(,(,+0.00 
■""-0.00 


iJ-uy.o 09 


0, 


(,(,+0.00 
•'-"J-0.00 


0, 


(,(,+0.00 
■'-"J-0.00 


0, 


(,(,+0.00 

■'J'-'-o.oo 


0, 


(,(,+0.00 
■'-"J-0.00 


0, 


(,(,+0.00 

■'J'-'-o.oo 



Table 13. Same as TablefTOl 



HD101412 HD104237 



SS73 44 



HD139614 



HD142666 HD142527 HD144432 











94.79 




73.75 




89.14 




26.95 




26.87 




701.54 




21.25 


Am. 01. 


0, 


1 fim 


0, 


(,(,+0.00 
■""-0.00 


0, 


(,(,+0.00 
"'^-0.00 


54, 


04+0.59 

■^^-0.54 


1 7 9c;+8.03 
-l-'.^0_14,91 


0, 


Og+1,85 
""-0,06 


0, 


(,(,+0.00 
■""-0.00 


0, 


nc:+0,98 
■""-0,05 


Am. 01. 


2, 


,0 fim 


0, 


(,(,+0.00 

■"^"-0.00 


0, 


(,(,+0.00 
■"'^-0.00 


0, 


(,(,+0.00 
■""-0.00 


47+14.20 
^■^' -7.65 


0, 


-.1+1,93 
■^^-0,10 


0, 


(,(,+0.00 
■""-0.00 


0, 


£.(5+2. 80 
■""-0.53 


Am. 01. 


5, 


,0 ^im 


0, 


(,(,+0.00 
■'^"-0.00 


0, 


(,(,+0.00 
■"'^-0.00 


0, 


(,(,+0.00 
■""-0.00 


0, 


-, 9+0.83 
■^^-0.11 


0, 


90+310 
■^"-0.20 


0, 


(,(,+0.00 
■""-0.00 


0, 


(,0+1.47 
■"'^-0.08 


Am. Py. 


0, 


,1 fim 


0, 


(,(,+0.00 
■'^"-0.00 


90, 


qc + 0.16 
■"'^-0.15 


43, 


K7+0.58 
■^' -0.56 


0, 


7Q+3.53 
■'^-0.77 


32, 


09+6.63 
■■^^-6.12 


93, 


90+0.23 

■^"-0.29 


95, 


40+0.60 


Am. Py. 


2, 


,0 fim 


88, 


-I 7+0.21 
■-'- ' -0.20 


0, 


(,(,+0.00 
■"'^-0.00 


0, 


(,(,+0.00 
■""-0.00 


69, 


00+2.01 

■■-"^-Les 


59, 


74+6.70 
■'^-6.98 


0, 


(,(,+0.00 
■""-0.00 


0, 


(,(,+0:00 
■""-0.00 


Am. Py. 


5, 


,0 fim 


0, 


ni +0,03 

■'J^-0.01 


0, 


(,(,+0.00 
■'-"^-0.00 


0, 


00+0.31 
■"^-o.os 


0, 


(,(,+0.15 
■""-0.00 


0, 


(,(,+0.16 
■""-0.00 


0, 


,(,+0.14 
■^"-0.09 


0, 


(,4+0.21 

■"^-0.03 


Fors. 


0, 


,1 /im 


0, 


00+0. 09 
■°'-'-0.08 


3, 


74+0.07 
■'^-0.06 


0, 


(,(,+0.00 
■""-0.00 


0, 


(,(,+0.00 
■""-0.00 


0, 


(50+0.13 

■"■'-0.14 


1, 


(,(,+0.09 
■'^"-0.09 


0, 


Qo+0.08 
■^"-0.09 


Fors. 


2, 


,0 fim 


2, 


nc+0.16 
■^"-0.13 


1, 


99+0.13 
■^^-0.11 


0, 


(,(,+0.03 
■""-0.00 


0, 


(,(,+0.01 
■""-0.00 


0, 


(59+0.28 
■"^-0.26 


4, 


70+0.23 
■ '"-0.22 


0, 


40+0.11 
■^"-0.13 


Enst. 


0, 


,1 fim 


3, 


9C+0.30 
■^"-0.32 


2, 


-,0+0.08 
■^"-0.08 


0, 


(,(,+0.00 
■""-0.00 


0, 


(,(,+0.00 
■""-0.00 


0, 


n'j+0.25 
■""-0.03 


0, 


(,(,+0.00 
■""-0.00 


0, 


-I 0+0.21 
■^"-0.11 


Enst 2, 


,0 fim 


2, 


74+0.39 
■ '^-0.37 


0, 


(,(,+0.00 
■"'^-0.00 


0, 


(,(,+0.00 
■""-0.00 


3, 


-,(,+0.42 
■^"-0.37 


5, 


90+0.31 
■^"-0.35 


0, 


(,(,+0.00 
■""-0.00 


1, 


7-, +0.20 
■ ' -'--0.28 


Silica 0, 


.1 fim 


2, 


70+0.08 
■ '"-0.09 


2, 


cn+0.05 

■""-0.04 


1, 


41+0.04 

■^-'--0.04 


0, 


00+0. 09 
■*»-0.08 


0, 


7c:+0.20 
■'"-0.21 


0, 


qn+0,04 

■""-0,04 


0, 


40+0.09 
■^'^-0.11 


Silica 2, 


.0 fim 


0, 


(,(,+0.00 
■'J'^-0.00 


0, 


(,(,+0.00 
■'-"J-0.00 


0, 


(,(,+0.00 
■""-0.00 


0, 


(,(,+0.00 
■""-0.00 


0, 


nQ+0.50 

■""-0.03 


0, 


(,(,+0,00 
■""-0,00 


0, 


(,-,+0.18 

"J^-o.oi 


Silica 5, 


,0 fim 


0, 


(,(,+0.22 

■'J"-o.oo 


0, 


(,(,+0.00 
■'-'"-0.00 


0, 


(,(,+0.00 
■""-0.00 


0, 


(,(,+0.36 
■""-0.00 


0, 


09+1.28 

■"^-0.30 


0, 


(,(,+0,00 
■""-0,00 


0, 


(,(,+0.00 
■""-0.00 
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Table 14. Same as TablefTOl 



HD144668 HD145263 HD150193 HD152404 51 Oph HD163296 VV Scr 











63.15 


30.21 




94.90 


6.79 




101.24 




31.49 




82.00 


Am. Ol. 


0. 


1 fim 


25, 


04+0,88 

'^^-0,88 


"■'-'"-0.00 


0, 


(,(,+0.00 
■"'-'-0.00 


18.66j:i74i 


0, 


(,(,+0.00 

■"'-'-0.00 


0, 


,001°: 


,00 
,00 


51, 


40+0.94 

■^^-0.98 


Am. Ol. 


2, 


fim 


0, 


„„+0.00 
■'^^-0.00 


no+°°° 
"■""-0.00 


0, 


(,(,+0.00 
■'J^J-o.oo 


r,,-, 01+20:10 
ZU.Oi_20.lO 


23, 


0.1+1.50 

■J •^-1.38 


0, 


,oo«: 


,01 
,00 


0, 


(,(,+0.00 
■'^'J-0.00 


Am. Ol. 


5, 


,0 /im 


0, 


(,(,+0.00 
■'J'^-0.00 


00+"°° 


0, 


(,(,+0.00 
■'^'-'-0.00 


'-'■'J^-0.00 


0, 


(,(,+0.04 

■'J^-o.oo 


0, 


■3itg: 


,26 
,31 





07+1.31 

•^'-0.27 


Am. Py. 


0, 


,1 fim 


70, 


cq+0.S9 
■°''-0.85 


o.oilHl 


68, 


f.o+0.96 


4R 4(,+ 10.58 
4D.4U_22.4g 


69, 


4(,+ 1.37 
■^'J-1.37 


64, 


,94t| 


,58 
,44 


31 


Qo+1.22 
■^''-1.08 


Am. Py. 


2, 


,0 fim 


0, 


(,(,+0.00 

■'J'-'-o.oo 


81.29l?.^?i 


18, 


70+0.91 
■ '•^-0.95 


6.80i^*7r 


0, 


(,(,+0.00 

■'J'-'-o.oo 


28, 


,86+3 


,53 
,13 


11, 


07+1.30 
■'J'-1.59 


Am. Py. 


5, 


.0 /im 


0, 


14+0.49 


10.99ti^of 


0, 


(,1+0.04 

■'J^-o.oi 


O.OOt'oTo 


0, 


(,(,+0.00 

■'^'J-o.oo 


0, 


,oo«: 


,00 
,00 


0, 


(,(,+0.00 
■'-"^-0.00 


Fors. 


0, 


,1 /im 


1, 


04+0.04 
■■^^-0.04 


1 74+0.21 


3, 


(,(,+0.05 
■""-0.05 


^■^"-0.16 


0, 


(,(,+0.00 
■'J'J-0.00 


1, 




12 

,12 





(,(,+0.04 

•^'J-o.oo 


Fors. 


2, 


,0 /im 


0, 


(,(,+0.00 

■^'J-o.oo 


^•"'5-0.40 


2, 


47+0.09 
-0.09 


n 61 +0-23 
'^■"^-0.21 


4, 


90+0.11 

■^^-0.11 


1, 


,96«: 


14 

,18 


0, 


00+0.08 
■°^-0.09 


Enst. 


0, 


,1 /im 


1, 


(,0+0.11 

■'^■^-0.13 


p. 00+0.43 
O.oo_0.40 


3, 


(.(,+0.07 
■""-0.07 


"■•^"-0.16 


1, 


1 1+0.23 
■^^^-0.23 


2, 


■671°: 


13 

,15 


2 


00+0.10 
•■^■^-0.10 


Enst 2, 


,0 /im 


0, 


,(,+0.19 


'^•'-'^-0.00 


0, 


(,(,+0.00 

■'J^J-o.oo 


*J-^'^-0.03 


1, 


c:c:+0.29 
■ "J '^-0.29 


0, 


,oo«: 


,00 
,00 


0, 


(,(,+0.00 
■""-0.00 


Silica 0, 


,1 /^m 


0, 


01 +0.06 
■^'-'--0.05 


ni+''-2'^ 
"■"^-0.01 


3, 


44+0.08 
■^^-0.07 


'-'■'J^-0.01 


0, 


on+0.07 
"^^-0.07 


0, 


,oo«: 


,03 
,00 


2 


00+0.07 
•U°-0.07 


Silica 2, 


.0 /im 


0, 


(,(,+0.00 
■'J'J-0.00 


no+°°° 
"■""-0.00 


0, 


(,(,+0.00 

■""-0.00 


'^■■'-'--0.22 


0, 


(,(,+0.00 

■^*J-o.oo 


0, 


,oo«: 


,00 
,00 


0, 


(,(,+0.00 

■^'J-o.oo 


Silica 5. 


,0 /^m 


0, 


(,(,+0.00 

■'J'-'-o.oo 


'^•'-"J-o.oo 


0, 


(,(,+0.00 
■'^'-'-0.00 


"■'J'-3.03 


0, 


(,(,+0.00 
■'J'J-0.00 


0, 


,001°: 


,00 
,00 





(,(,+0.00 
•'-"J-0.00 



Table 15. Same as TablefTOl 



T CrA HD179218 WW Vul HD190073 HD203024 











15.84 




651.54 




69.70 




83.84 




197.33 


Am. 01. 


0, 


1 fim 


0, 


(,(,+0.01 
■'-"^-0.00 


0, 


(,(,+0.00 
■'^'-'-0.00 


0, 


(,(,+0.00 
""-0.00 


0, 


00+°°° 
""-0.00 


0, 


00+°°° 
""-0.00 


Am. 01. 


2, 


fim 


0, 


(,(,+0.00 
■'-"^-0.00 


0, 


(,(,+0.00 
■'^'-'-0.00 


0, 


(,(,+0.00 
■""-0.00 


0, 


00+°°° 
■""-0.00 


0, 


00+°-°° 
■""-0.00 


Am. 01. 


5, 


,0 ^im 


0, 


(,(,+0.00 
•^J'^-o.oo 


20, 


00+0. 89 
■"'-'-1.33 


0, 


(,(,+0.00 
■""-0.00 


0, 


00+°°° 
■""-0.00 


0, 


00+°°° 
■""-0.00 


Am. Py. 


0, 


,1 fim 


35, 


71 +2.36 
• ' ^-2.09 


66, 


74+1.26 
■ '^-0.91 


9, 


(,p.+3.88 
■""-'-3.73 


93. 


7C+0.23 
'°-0.18 


86. 


10+1.15 

■^*'-1.29 


Am. Py. 


2, 


,0 fim 


48 
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Fig. 1. — Comparison of scatterin g theories and differe n t sets of optical c onstants for crystalline 
forsterite (jServoin &: Pirioul ()l973l ): ISogawa et alJ (j2006l ): ISuto et alJ (j2006l )). The applied scatter- 
ing theories were, (a) Mie theory, (b) CDE, and (c) DHS with a maximum volume filli ng factor of 1.0 



For c omparison the optical efficiencies of forsterite measured on free-flying particles ([Tamanai et al 



2003 ) are shown. 




Fig. 2. — Same as Figuredl but for longer wavelengths. 
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Fig. 3. — Absorption efficiencies of the amorphous sihcates. Sohd, dotted and dashed hnes show 
the absorption efficiencies calculated from the optical constants (see Table[2]) using the DHS theory 
with a maximum volume ffiling factor of 1.0, 0.7 and 0.0, respectively. The zero ffiling factor 
solution is identical to the solution of Mie theory. 
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Fig. 4. — PAH emission profiles used for the analysis. All band profiles are normalized to their 
peak-value. For the details of the derivation of the different profiles and the notation of the profile 
names (X1-X6) see Sectioning] 
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Fig. 5. — Fits to the short wavelength range (5-17 /xm). The observed IRS spectrum is shown 
with black dots, while the red line shows the best-fit model. For each fit the residuals ((F™°dei_ 
pobs^/pobs X -^QQ-^ g^^g j^igQ shown. 
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Fig. 6. — Same as FigureO 
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Fig. 7. — Same as FigureO 
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Fig. 8.— 



Same as FigureEl 
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Fig. 9. — Fits to the long wavelength range (17-35 /xm). The observed IRS spectrum is shown 
with black dots, while the red line shows the best fit model. For each fit the residuals ((F™°dei_ 
pobs^/pobs X -^QQ-^ g^^g j^igQ shown. 
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Fig. 10. — Same as Figure^ 
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Fig. 11. — Same as FigureEl In the spectrum of HD139614 the emission feature between 26 and 
31 fim is apparent and caused by problems in the 15th and 12th order of the long-high module (see 
Section[22]). 
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Fig. 12. — Same as FigureEl 
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Fig. 13. — Left: grain shape effects in tlie case of a 0.1 /xm sized forsterite grain in the 10 /xm region. 

Right: the same as for amorphous sihcates with ohvine stoichiometry and with Fe / (Mg+Fe) = 0.5. 
It can be seen that differences in the calculated absorption efficiencies by DHS and CDE theories are 
far smaller, than between Mie-theory and the other two scattering theories. For crystalline silicates 
the differences between the calculated optical efficiencies by arbitrary two scattering theories are 
larger than a few percent, which is a typical error level in our fits. 



100 - . * 

* 

X - , • 



♦ • • r • • 



10 



0.01 0.10 



cryst 

Fig. 14. — Reduced of the fits as a function of crystallinity for the 7-17 /xm region. In general 
spectra with higher crystallinity have higher reduced in the fits suggesting possible weaknesses 
in our crystalline dust model. 
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Fig. 15. — Ratio of the normalized fluxes at 11.3 fim to that at 9.8 vs. peak- 
to-continuum ratio of the 10 fim sihcate comp lex. Normalized flux was calculated as 



pnorm 
1/ 



1 + (F' 



obs 
u 



pcont^^ < F^°°* >, according to Ivan Boekel et al.l (j2005l ). In the box (dashed 
lines) sources have the most pristine 10 /um silicate feature similar to that in the ISM. These sources 
were selected to test the amorphous dust population. 




Fig. 16. — Comparison of different datasets of amorphous silicates used for the fit. The three 
panels show the results for a) HD36112, b) HD152303 and c) HD144432, respectively. AMIXl and 
AMIX3 mixtures consist of iron- free silicates while for AMIX2 Fe/(Mg+Fe)=0.5. For the details 
of the different amorphous silicate mixtures, see Section HT2l The AMIXl gives always a lower 
than either of the two other mixtures. 
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Fig. 17. — Classification of the sources on the basis of the SED (jvan Boekel et al.ll2005l ). The 
plotted quantities are the ratio of the near-infrared to infrared luminosities vs. IRAS 12 fim 
color (mi2 - mgo = —2.5 x logFi2/-F6o)- The dashed line marks the boundary between Group / and 
Group //sources, according to Fnir/Fir = 1.5x(mi2 - meo)- 



t T 
• • 



□ Group la 
• Group I la 



10 12 



Fig. 18. — Mass-averaged grain size vs. disk flaring. The disk flaring is empiricany parameterized 
by the ratio of the flux densities at 24 fim and 8 /um. A trend is clearly visible within Group Ila, 
such that sources with steeper mid-infrared SED slope have larger grains in the disk atmosphere. 
In the case of the outliers in Group la the calculated flux ratio is not likely to measure the disk 
flaring only, but it is also influenced by other disk parameters (see the text for the details). 
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Fig. 19. — Abundance ratios of enstatite and forsterite derived from the short and long wavelength 
flts for (a) Group I and (b) for Group II sources. The dashed line marks the 1:1 ratio between 
the two plotted quantities. It can be seen that forsterite is more abundant than enstatite at longer 
wavelengths (i.e. outer disk) while at shorter wavelengths (i.e. inner disk) the situation is the 
opposite. 
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Fig. 20. — Com parison of the MACs of enstatite derived from the Spitzer IRS spectra and laboratory 
measurements (jChihara et al.ll2002l ). See Section lS^l for the details of the derivation of the enstatite 
absorption coefficients from the IRS spectra. It can be seen that the absorption coefficients derived 
from the IRS spectra are the most similar to that of "En90", which is a crystalline pyroxene with 
about 10 % iron content. 
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Fig. 21. — Comparison of the MACs of forsterite derived from the Spitzer IRS spectra (see Sec- 
tion l5.3p . calculated using DHS theory (used for spectral decomposition) and meas ured in labora- 



tory. In the laboratory experiment the MACs were measured on free flying particles (jTamanai et al 



20091 ). Although the positions of the bands at 10 //m and 11.3 //m are better reproduced by DHS 
calculations than laboratory measurements, in the case of the 16 //m band the situation is the 
opposite. Neither of the two MAC curve (DHS, laboratory measurement) can reproduce all the 
observed peak positions of forsterite in the same time. 



